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I. 


INTRODUCTION 


The  objective  of  this  study  is  to  determine  the  relative  importance 
of  carburization  to  the  overall  degradation  of  a  gun  barrel  by  propellant 
gases.  Carbon  monoxide,  the  pre-eminent  carburizing  atmosphere,  is  present 
in  all  existing  gun  propellants,  in  quantities  large  enough  to  suggest  a 
significant  contribution  to  the  thermo chemical  history  of  a  ballistic  system. 
Carburization  is  normally  thought  of  as  a  steady-state  diffusion  process  which 
enhances  a  steel's  compressive  strength,  surface  hardness  and  wear  resistance. 
But  as  a  possible  internal  ballistic  process,  carburization  can  only  be  present 
as  a  transient,  on-the-surface  phenomenon,  since  high  propellant  gas  tempera¬ 
tures  cannot  increase  the  rate  of  carbon  diffusion  to  a  level  where  subsurface 
microstructure  is  noticeably  altered  in  the  span  of  several  milliseconds. 
Furthermore,  there  is  no  regenerative  medium  present  in  a  ballistic  system  to 
dissociate  the  CO2  produced  in  the  Fe-CO  reaction.  Any  CO2  produced  by  CO 
induced  Fe3C  formation  is  either  expelled  by  forced  convection  or,  more 
importantly,  reverses  the  carburizing  process  as  barrel  temperature  decreases. 
It  is  this  cyclical  history  of  carburizing  and  decarburizing  of  the  gun  tube 
surface  during  a  ballistic  event  that  is  believed  to  be  responsible  for  any 
erosion  potential  exhibited  by  the  presence  of  carbon  monoxide.  It  is  further 
postulated  that  in  high  performance  guns,  where  barrel  surface  temperatures 
already  hover  near  the  steel !s  normal  solidus  temperature,  an  increase  in  car¬ 
bon  concentration,  even  for  a  few  milliseconds,  will  artificially  lower  the 
melting  point  to  where  the  barrel fs  surface  shear  strength  is  lowered  or  lost. 
This  program  has  been  designed  to  ascertain  the  significance  of  these  pro¬ 
cesses  and  quantify  their  contribution  to  gun  barrel  erosion  and  cracking. 

Research  is  being  performed  by  utilizing  Calspan's  Shock  Tube  Gun 
(STG)  facility,  which  duplicates  the  thermodynamics  of  gun  chamber  combustion 
with  a  tube  chamber  incorporating  a  moving  wall,  i.e.,  the  leading  face  of  a 
gas  driven  piston.  The  gases  compressed  in  the  tube  chamber  by  piston  motion 
vent  in  a  conventional  manner  by  forcing  a  projectile  down  the  adjoining  gun 
barrel.  During  the  venting  cycle,  the  gases  must  pass  through  the  nozzled 
flow  channel  of  a  gun-steel  sample  positioned  in  front  of  the  barrel  entrance. 
As  a  result,  the  channel  wall  of  the  sample  experiences  histories  of  pressure, 
temperature  and  forced  convective  heating,  similar  to  those  of  a  gun  barrel. 
Unlike  the  combustion  system  where  propellant  formulation  dictates  a  set 
flame  temperature  and  impetus,  the  STG  permits  independent  variation  of  any 
or  all  conditions  affecting  ballistic  histories. 

Several  computer  programs  provide  analytic  support  for  the  STG  tests. 
One  code  models  the  STG  gas  compression  cycle  using  a  van  der  Waals  equation  of 
state.  It  computes  gas  pressure  and  temperature,  convective  heating,  total 
heat  input  and  in-wall  temperatures  of  the  test  sample  surface.  This  program 
is  used  mainly  in  a  predictive  mode  to  select  gas  composition  and  peak  pres¬ 
sure  for  an  actual  test.  Another  code  accepts  gas  mixture  and  pressure  data 
from  tests  conducted  and  computes  equilibrium  temperature  and  species  concen¬ 
tration  values  along  the  correct  isentropes.  This  information  can  then  be 
returned  to  the  preceding  program  to  improve  its  prediction  of  barrel  heating 
and  wear.  STG  testing  and  computer  analysis,  as  performed  in  this  program,  are 
intended  to  expand  present  knowledge  of  wear  and  erosion  phenomenology. 
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II.  SHOCK  TUBE  GUN  FACILITY 


2.1  STG  Concept 

The  Shock  Tube  G_un>  as  designed  and  developed  by  Calspan  applies 
shock  tube  principles  to  the  study  of  interior  ballistics.  The  facility 
consists  of  a  driver  gas  chamber,  a  driven  tube  containing  a  latchable 
flying  piston,  an  instrumented  gas  collection  chamber  and  an  instrumented 
gun  tube  containing  a  projectile. 

One  unique  design  feature  of  the  STG  is  the  driven  piston  which, 
by  its  presence,  affords  physical  isolation  of  the  driver  gas,  normally 
nitrogen,  and  the  chamber  or  test  gas,  which  varies  in  composition  with  test 
objectives.  The  driven  piston,  by  virtue  of  its  mass  in  conjunction  with 
the  projectile’s  mass,  also  controls  the  compression  history  of  the  test 
gas,  enabling  the  facility  to  duplicate  the  interior  ballistic  environment 
of  various  large  caliber  guns,  up  to  and  including  an  8”  howitzer. 

Another  key  design  feature  of  the  STG  is  the  test  section  contained 
within  the  gas  collection  chamber.  It  accepts  specimens  of  various  physical 
configurations  equipped  with  in-wall  thermocouples  and/or  surface  heat  flux 
meters,  and  holds  them  in  place  adjacent  to  the  attached  barrel.  In  effect, 
this  design  offers  a  replaceable,  highly  instrumented  bore  entrance  to  an 
otherwise  conventional  gun  tube. 

With  suitable  variation  in  test  parameters,  made  possible  by  the 
design  of  the  STG,  including  driver  pressure,  piston  mass,  test  gas  composi¬ 
tion,  specimen  configuration,  shot  start  capability  and  projectile  mass,  one 
may  investigate  any  or  all  of  the  internal  factors  affecting  ballistic  pheno¬ 
mena,  such  as  the  gun  barrel  erosion  and  cracking,  which  the  present  study 
addresses . 


2.2  Construction 


Table  1  lists  the  present  structural  dimensions  of  the  STG  as  used 
in  this  study,  and  dictated  by  adiabatic  compression  modeling  and  material 
availability. 

As  shown  in  Figure  1,  the  projectile  launch  components  consist  of 
the  191  millimeter  driven  tube,  the  instrumented,  high  pressure,  test  gas 
collection  chamber  and  a  30mm  smooth-bore  barrel.  These  are  supported  on  a 
shock  table  which  is  free  to  move  on  tracks  in  the  direction  of  piston 
motion,  during  the  severe  impulse  loading  caused  by  unbalanced  -chamber  pres¬ 
sure  due  to  test  gas  compression  by  the  decelerating  piston.  This  floating 
mount  system  minimizes  shearing  forces  to  the  supporting  base  structure, 
but  requires  an  adjustable  pneumatic  brake  on  the  driven  tube  to  absorb 
the  axial  loading  on  the  launch  components,  primarily  the  driven  tube  itself. 
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Table  1. 


Shock  Tube  Gun  Characteristics 


Configuration  Data: 


Driven  Tube  I.D. 

0.191m 

(7.5  in.) 

Driven  Tube  Length 

24.6m 

(970  in.) 

Piston  Area 

0.0285m2 

(44.179  sq. 

in.) 

Piston  Mass 

Up  to  91  kg. 

(200  lb.) 

Projectile  Diameter 

30mm 

(1.181  in.) 

Projectile  Area 

706mm^ 

(1.095  sq. 

in.) 

Projectile  Mass 

Up  to  .91  kg. 

(2  lb.) 

Driver  Volume 

0.885m3 

(54,000  cu. 

in.) 

Chamber  Volume 

2140mm3 

(130.8  cu. 

in.) 

Pressure  -  at  release  of  projectile 

Barrel  Length 

Variable 

4.57m 

(180  in.) 

The  projectile  capture  components  consist  of  a  telescoping  tube 
coupled  to  the  barrel  muzzle,  a  projectile  blast  chamber,  and  a  sand  filled 
tube  to  decelerate  and  catch  the  projectile.  The  telescoping  tube,  the 
purpose  of  which  is  to  permit  independent  motion  of  the  shock  table  and 
blast  chamber,  contains  replaceable  screens  for  measuring  projectile  velocity. 
The  blast  chamber  reduces  the  noise  and  pressure  levels  as  the  projectile 
exits  the  barrel . 

Figure  2  shows  the  chamber  and  toggle  restraint  system  needed  to 
contain  the  high  chamber  pressures  and  associated  axial  loads.  Chamber  pres¬ 
sures  are  sensed  using  piezoelectric  transducers.  The  entrance  region  of 
the  launch  tube  can  accommodate  pressure,  heat  flux,  and  erosion  sensing 
devices . 


The  piston,  which  is  used  to  compress  the  test  gas,  is  made  from 
4340  steel  and  weighs  68  kg.  including  the  latching  block  on  its  rear  face 
which  secures  it  at  the  upstream  end  of  the  driven  tube,  prior  to  release 
or  "firing."  Gas  seal  is  maintained  using  "T"  rings  at  the  front  and  rear  of 
the  piston.  Three  brass  wear  rings  or  "bore  riders"  are  used  to  prevent 
steel-to-steel  contact  between  piston  and  tube.  A  buffer  projection  on  the 
face  of  the  piston  and  a  complementary  piston  stop  ring  at  the  downstream 
end  of  the  driven  tube  prevent  direct  impact  of  the  piston  into  the  test 
gas  collection  chamber  in  the  event  that  the  compressed  test  gas  develops 
insufficient  pressure  or  loses  pressure  prematurely  due  to  chamber  seal 
failure . 
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191  mm  DRIVEN 
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Figure  1  PROJECTILE  LAUNCH  AND  CAPTURE  COMPONENTS  OF  THE  SHOCK  TUBE  GUN 
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Figure  2  PISTON  CHAMBER  AND  TOGGLE  RESTRAINT  COMPONENTS 
OF  THE  SHOCK  TUBE  GUN 
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III.  TEST  PREPARATIONS 


3.1  Specimen  Selection 


The  primary  objective  of  this  study  was  to  determine  if  certain 
propellant  gas  conditions  enhance  barrel  erosion  and  cracking.  To  correlate 
the  test  data,  the  specimens  used  were  made  nearly  identical  in  shape  and 
composition,  i.e.,  a  4340  steel  cylinder,  38.1mm  in  length,  31.75mm  in 
diameter,  and  bored  concentrically  to  12.7mm,  to  create  a  sonic  flow  condi¬ 
tion  during  the  tests.  The  flow  channel  inlet  was  radiused  to  reduce 
turbulence  and  to  increase  heat  flux  over  a  larger  portion  of  the  flow 
channel  surface.  The  samples  were  small  enough  to  fit  in  the  specimen 
stage  of  Calspanfs  Scanning  Electron  Microscope  (SEM)  such  that  bore  sur¬ 
face  examinations  could  be  conducted  without  using  replicas.  Conversely, 
the  samples  were  large  enough  to  register  finite  changes  in  weight  and 
bore  diameter,  resulting  from  test  conditions.  Mass  changes  were  measured 
in  tenths  of  a  milligram  on  an  analytical  balance.  Diametral  recession  was 
measured  to  within  10~4mm  at  four  specific  axial  locations,  as  shown  in 
Figure  3.  Also  shown  in  the  figure  are  the  ports  for  in-wall  thermocouples 
which  are  used  to  determine  the  integrated  heat  input. 


4.76R 
- X 


38.1 


Figure  3  SHOCK  TUBE  GUN  TEST  SAMPLE 
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3.2 


Heat  Transfer  Instrumentation 


A  primary  measurement  of  the  study  is  the  amount  of  bore  heating 
associated  with  each  test.  For  this  measurement,  two  in-wall  thermocouples 
were  installed  in  each  sample,  at  distances  approximately  0.5mm  from  the  bore 
surface.  The  method  of  installation  is  shown  in  Figure  4.  Each  of  these 
thermocouples  independently  may  be  used  to  determine  net  heating  to  the  bore. 
Total  heat  input  is  calculated  from  the  in-wall  thermocouple’s  output  by 
use  of  methods  developed  and  reported  by  Calspan.1 2  Briefly,  conversion  of 
thermocouple  output  (millivolts  vs.  time)  to  total  net  heat  input  per  unit 
area  is  made  by  use  of  the  relation 


Q(t)  =  AT  (t)  l/rrkcpt  (1) 

where  Q(t)  is  the  net  bore  heat  input 

AT(t)  is  the  indicated  change  in  in-wall  temperature  as  a  function 
of  time 

k  is  the  thermal  conductivity 
cp  is  the  heat  capacity  per  unit  volume 
t  is  the  time  after  start  of  heating. 

Data  reduction  procedure  consists  of  calculating  Q(t)  using 
Equation  (1)  at  successive  time  intervals  to  produce  a  curve  of  Q(t)  vs.  t 
which  becomes  asymptotic  to  the  true  heat  input.  To  compute  the  correct 
asymptote,  time  zero  for  the  start  of  heating  must  be  established  accurately 
on  the  thermocouple  trace.  It  was  concluded  in  the  previous  STG  study2>  that 
correctly  shaped  asymptotes  were  most  consistently  produced  by  placing  time 
zero  at  the  intersection  of  the  trace  baseline  and  the  average  slope  of  the 
initial  heat  pulse’s  leading  edge. 

The  Table  2  values  for  Heat  Input  are  taken  from  end  points  on 
the  Q(t)  vs.  t  curves  that  have  been  divided  by  a  correction  factor  to 
account  for  the  heat  flux  being  dispersed  over  an  increasing  surface  area 
as  it  passes  through  the  test  specimen’s  radial  wall.  This  correction  factor, 
p(r),  based  on  geometric  considerations,1  is  calculated  from  the  relation 

n(r)  =  1  -0.32e'3'71r 


For  the  test  specimen  bore  radius  of  6.35mm,  p(r)  =  .87 


1.  F.A.  Vassallo,  ’’Mathematical  Models  and  Computer  Routines  Used  in 
Evaluation  of  Caseless  Ammunition  Heat  Transfer,”  Calspan  Report  No. 
GM-2948-Z-1,  June  1971. 

2.  E.B.  Fisher,  C.C.  Morphy,  ’’The  Role  of  Oxygen  in  Gun  Barrel  Erosion  and 

Cracking--A  Shock  Tube  Gun  Investigation,”  U.S.  Army  Ballistic  Research 
Laboratory,  ARRADCOM,  Contract  Report  No.  ARBRL-CR-00427,  April  1980. 
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SHOCK  TUBE  GUN  TEST  SECTION 


TEST  SPECIMEN 


SILICONE 

GREASE 


THERMOCOUPLE 
WIRES  (40  GAGE) 
CHROMEL  ALUMEL 

w 


BORE  SURFACE 


Figure  4  IN-WALL  THERMOCOUPLE  INSTALLATION 
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3.3 


Selection  of  Test  Conditions 


The  program  objective  was  to  experimentally  characterize  the  role 
of  carburization  in  gun  barrel  erosion.  To  insure  meaningful  results,  the 
selection  of  test  conditions  had  to  meet  three  basic  criteria.  First,  all 
but  the  "inert"  test  gas  mixtures  contained  carbon  monoxide  to  duplicate 
the  reducing  species  in  propellant  gas  atmospheres,  particularly  those  with 
low  flame  temperatures,  such  as  nitramines  which  have  relatively  high  CO 
concentrations.  Test  concentrations  of  carbon  monoxide  were  on  the  average 
higher  than  any  found  in  actual  propellant  gas,  but  were  deemed  appropriate 
for  this  initial  study  of  the  gross  effects  of  carburization  in  relation  to 
other  ballistic  phenomena  such  as  oxidation  or  pure  melting  which  were  pre¬ 
viously  shown  to  influence  barrel  erosion  and  cracking. 2,3  The  logical  exten¬ 
sion  of  this  present  study  would  be  an  investigation  of  quantitative  propellant 
gas  chemistry,  incorporating  all  of  the  phenomena  reported  on  to  'date. 

The  second  criterion  to  be  met  was  the  selection  of  thermodynamic 
conditions  indigenous  to  large  caliber  guns.  To  a  large  extent  the  interior 
conditions  were  dictated  by  the  STG's  present  physical  configuration,  including 
the  mass  of  the  driven  piston  which  transfers  energy  from  the  compressed  gas 
to  the  test  mixture.  However,  the  mass  of  the  test  mixture  often  varied 
from  run  to  run  and  appropriate  calculations  were  made  with  the  use  of  the 
STG  computer  code  to  determine  the  correct  driver  pressure  that  would  produce 
the  desired  real  gun  conditions,  namely  pressures  in  the  neighborhood  of 
300  MPa  and  temperatures  near  the  melting  point  of  the  test  specimen  alloy, 

4340  steel. 

A  third  criterion  for  the  selection  of  test  conditions  was  concerned 
with  event  time  variation.  Since  carburization  is  dependent  upon  time  in 
addition  to  concentration  gradient  and  temperature,  an  effort  was  made  to 
alter  the  profile  of  the  heat  and  pressure  pulse  while  maintaining  its 
peak  amplitude  at  a  constant  value  through  several  tests.  It  was  felt  that 
if  carburization  affected  gun  barrel  wear  within  a  period  of  several  milli¬ 
seconds,  then  a  variation  in  time  of  several  fractions  of  a  millisecond  would 
also  effect  barrel  wear. 

In  meeting  the  above  criteria,  the  test  conditions  selected  were 
able  to  separate  the  contributions  of  chemistry,  pressure,  temperature  and 
time  to  carburization  potential  and,  in  turn,  its  effect  upon  barrel  erosion. 


3.  F.A.  Vassallo,  W.R.  Brown,  "Shock  Tube  Gun  Melting  Erosion  Study,"  USA 
Ballistic  Research  Laboratory,  ARRADCOM,  Contract  Report  ARBRL-CR-00406 , 
January  1979.  (AD  #A076219) 
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IV.  TEST  PROCEDURE  AND  RESULTS 


4.1  Test  Procedures 


Operation  of  the  STG  to  collect  test  data  regarding  ballistics, 
heating  and  erosion  followed  a  fixed  experimental  pattern.  Prior  to  a  run 
being  conducted,  the  facility  was  inspected  for  damage  from  the  preceding  test 
and  needed  repairs  were  performed.  Components  such  as  the  driven  tube,  piston, 
and  gas  collection  chamber  which  contact  the  test  gas  during  a  run  were  care¬ 
fully  cleaned  to  eliminate  contaminants  from  their  surfaces.  All  wearing 
surfaces,  including  piston  "0"  rings  and  bore  riders,  tube  seals  and  stop 
ring  buffers  were  replaced  if  their  wear  limits  were  reached.  The  piston  was 
then  inserted  in  the  upstream  end  of  the  tube  and  latched  to  the  driver 
release  mechanism.  A  new  projectile  was  inserted  in  the  barrel. 

The  numbered  specimen  was  weighed  prior  to  testing  using  an  analyti¬ 
cal  balance  for  initial  mass,  given  a  final  cleaning  with  freon,  and  then 
installed  in  the  sample  holder  within  the  gas  collection  chamber.  M-ll  mechani 
cal  pressure  gauges  were  also  installed  in  the  chamber.  Thermocouples  were 
inserted  through  the  chamber  wall,  positioned  in  the  sample  wall  and  then  the 
chamber/barrel  assembly  was  lowered  into  position  and  sealed  to  the  downstream 
end  of  the  driven  tube  with  hydraulic  toggles. 

After  installation  of  the  projectile  and  specimen,  the  entire  tube/ 
chamber  cavity  was  evacuated  to  a  pressure  of  2.0mm  Hg  or  less.  If  vacuum 
was  maintained  for  a  reasonable  time,  indicating  seal  integrity,  the  cavity 
was  purged  with  argon,  re-evacuated,  again  purged  with  argon,  and  evacuated 
for  a  third  time.  The  cavity  was  then  charged  to  the  local  atmospheric 
condition  with  the  required  partial  pressures  of  the  gases  selected  for  the 
test  mixture.  These  partial  pressures  are  dependent  upon  the  mix  ratio  desired 

Equations  for  establishing  partial  pressure  settings  were  derived 
from  the  Dalton  model  of  ideal  gas  mixtures,  which  assumes  the  following: 

1.  The  moles  of  mixture,  n,  equals  the  sum  of  the  moles 
of  the  component  gases,  nA  +  nN  +  nco,  where  A,  N  and 
CO  are  subscripts  referring  to  argon,  nitrogen  and 
carbon  monoxide,  respectively. 

2.  Each  component  gas  in  the  mixture  occupies  the  entire 
mixture  volume,  V,  which  in  this  case  is  the  'volume 
of  the  driven  tube. 

3.  The  temperature,  T,  of  the  components  before  and 
after  mixing  remains  constant. 

4.  The  mixture  pressure,  P,  in  this  case,  1  atmosphere, 
is  reasonably  low,  to  assure  near  ideal  gas  behavior. 
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For  the  components: 


For  the  mixture: 


PAV  =nART 
PNV  =  nNRT 

pcov  =  ncoRT 


PV 


nRT 


Since  V/RT  is  a  constant  in  all  the  equations: 


P  P  P 

__A  _  __N  JX>  =  P 

nA  nN  nC0  n 


Rewriting: 


PA 

nm 

A  _ 

A 

P 

n 

N 

_N 

P 

n 

pco 

 nco 

p 

n 

That  is,  for  each  component  of  a  mixture  of  ideal  gases,  the  mole  fraction 
and  the  ratio  of  the  partial  pressure  to  the  total  pressure  are  equal. 

Upon  completion  of  test  gas  charging,  the  mixture  was  given  time 
to  equilibrate  in  the  driven  tube  while  the  required  instrumentation  including 
pressure  transducers,  thermocouples  and  velocity  screens  were  connected  to 
suitable  recording  devices  and  checked  for  correct  operation.  The  piezoelec¬ 
tric  pressure  output  was  stored  for  future  playback  at  a  tape  speed  of  .76  mps 
(30  ipsj  on  a  Bell  §  Howell  7  channel  recorder.  The  thermocouple  output  was 
recorded  on  a  CEC  oscillograph  at  a  paper  speed  of  5  ips .  The  projectile  tran¬ 
sit  time  between  velocity  screens  was  measured  directly  on  a  Tektronix  oscillo¬ 
scope  with  trace  storage  at  1  ms  per  division. 

If  the  instrumentation  checked  out  satisfactorily*,  the  nitrogen  tank 
farm  valve  was  opened,  the  tube  air  brake  was  charged,  the  driver  was  pres¬ 
surized  to  the  desired  level  for  the  experiment,  recording  devices  were 
activated  and  the  piston  was  released. 

After  exhausting  residual  driver  pressure,  the  air  brake  was  bled, 
and  the  chamber/barrel  assembly  was  decoupled  from  the  tube.  The  specimen 
was  carefully  removed,  inspected,  weighed  and  measured  diametrally  at  the  afore¬ 
mentioned  axial  locations.  Hard  copy  was  made  of  all  test  data  for  further 
reduction  and  analysis. 
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4.2 


Test  Matrix 


A  total  of  32  tests  were  conducted  during  this  program.  A  summary 
of  the  test  conditions,  mass  loss,  the  recession  measurements  along  the  length 
of  the  sample’s  flow  channel  and  total  heat  input  as  measured  by  the  in-wall 
thermocouples  and  corrected  for  radial  dispersion,  are  presented  in  Table  2. 

The  test  matrix  was  divided  into  essentially  five  groups.  The  first 
group,  tests  63  through  68  established  a  wear  curve  for  a  10  percent  carbon 
monoxide  mixture,  at  a  fixed  driver  pressure.  The  remaining  test  gas  cons¬ 
tituents,  nitrogen  and  argon  were  varied  to  alter  the  ratio  of  specific 
heats,  which  was  the  sole  independent  variable  altering  flow  conditions  between 
tests.  Two  subsequent  tests,  69  and  70  adding  5  percent  hydrogen  to  the  first 
group’s  formula,  were  conducted  to  ascertain  the  peak  flow  conditions  obtained 
with  a  ’’lighter”  gas  mixture  with  the  same  concentration  of  active  constituents 
as  tests  63  and  64.  Test  71  was  conducted  as  an  ’’inert”  comparison  of  tests  63 
and  69,  since  the  peak  gas  temperatures  experienced  with  the  first  group  were 
well  above  threshold  of  erosion  values,  in  a  temperature  region  where  no 
previous  STG  data  had  been  collected  and  where  no  baseline  melting  had  been 
quantified. 

The  second  group  of  runs,  tests  72  through  75  and  also  77,  better 
defined  the  inert  wear  curve,  begun  in  the  previous  program^,  near  and  above 
the  erosion  threshold  for  flow  conditions  experienced  in  existing  large 
caliber  guns,  using  hot  double-base  propellants. 

The  third  group  of  runs,  tests  78  through  82,  established  a  wear 
curve  for  high  CO  concentration  mixtures,  by  direct  substitution  of  45.5 
percent  carbon  monoxide  for  an  equal  mole  fraction  of  nitrogen  used  in  the 
previous  group  of  ’’inert”  tests.  Since  CO  and  N2  have  the  same  molecular 
weight,  peak  conditions  of  tests  performed  at  identical  driver  pressures 
varied  little  between  the  second  and  third  groups.  Direct  comparison  of 
test  results  of  these  two  groups  proved  of  primary  importance  in  establishing 
the  enhanced  erosivity  of  propellant  gases  containing  large  quantities  of  CO. 

A  fourth  group  of  runs,  tests  83  through  88,  attempted  to  establish 
any  concentration  threshold  for  CO  that  may  exist  between  the  ’’inert”  no  CO 
atmosphere  of  group  2  mixtures  and  the  high  CO  concentration  mixtures  of  group 
3.  The  combined  fractional  percentage  of  CO  and  N2  remained  constant  at  45.5 
The  mass  of  the  mixture,  as  previously  explained,  also  remained  constant  for 
all  tests  in  this  group. 

A  fifth  and  final  group  of  runs,  tests  89  through  94,  were  conducted 
for  purposes  of  comparison  with  the  results  of  group  3  (45.5%  CO-no  N2)  tests. 
4%  He  was  substituted  for  an  equal  mole  fraction  of  argon,  which  ’’lightened” 
the  gas  mixture,  but  left  its  chemical  ’’activity”  unchanged.  The  purpose  of 
'these  tests  was  to  ascertain  the  contribution  of  heating  time  to  the  CO 
erosion  mechanism,  in  terms  of  rise  time  and  period  of  duration  of  the  heat 
pulse  above  a  critical  bore  temperature,  i.e.,  the  solidus  temperature  of  4340 
steel  with  appropriate  modification  for  local  (bore  surface)  chemistry  changes 
determined  from  comparison  of  group  3  results  and  STG  code  computations. 
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4.3 


Test  Gas  Mixture  Isentropes 


Adiabatic  compression  of  the  test  gas  in  the  STG  generates  the  pressure 
and  temperature  profiles.  The  test  gas  composition  alters  the  relationship 
between  the  gas  temperature  and  pressure  for  a  given  compression  ratio.  This 
relationship  is  shown  in  Figure  5  for  the  major  gas  compositions  used  in  this 
program.  Test  numbers  indicate  the  peak  position  of  each  run  on  its  respec¬ 
tive  isentrope.  The  isentropes  were  computed  by  an  equilibrium  combustion 
program.  Appendix  II  contains  a  brief  description  of  the  code  and  a  sample 
printout.  This  computer  program  was  also  used  to  determine  the  equilibrium 
composition  of  the  gas  mixture  at  peak  flow  conditions.  In  this  way,  with 
a  specified  gas  mixture  and  measured  pressure,  the  temperature  and  gas  cons¬ 
tituents  generated  as  a  result  of  the  dissociation  and  chemical  reactions 
were  determined  for  the  compression  cycle. 

It  is  important  to  note  that  the  equilibrium  code  treats  a  given  gas 
mixture  as  a  closed  system  with  no  chemical  or  thermal  interaction  with 
its  boundary.  The  code  does  not  recognize  the  presence  of  a  gun  barrel  wall 
that  is  removing  heat  from,  and  exchanging  molecules,  including  C,  with  the 
compressed  gas  mixture.  The  significance  of  gas-boundary  interaction  is 
demonstrated  by  the  isentropes  of  Figure  5.  The  45.5%  N2  (inert)  mixture 
reaches  the  solidus  temperature  of-  4340  steel  at  a  lower  pressure  than  the 
45.5%  CO  mixture,  indicating  a  lower  expected  threshold  of  erosion  for 
nitrogen.  As  will  be  shown  in  Section  V,  just  the  opposite  effect,  resulted 
in  actual  STG  testing. 

Table  3  lists  the  active  constituents  by  mole  fraction  of  each  test 
gas  mixture,  both  at  ambient,  pre-test  conditions  and  at  the  peak  pressure 
and  temperature  conditions  of  each  test.  Nitrogen,  though  classified  as 
"inert"  when  not  in  the  presence  of  carbon  monoxide,  has  been  included  in 
the  table  to  illustrate  its  potential  for  forming  oxides,  acids  and  cyanides 
with  active  mixture  components  during  a  test.  These  nitrogen  compounds,  all 
gaseous  byproducts,  contribute  little  to  direct  barrel  wear  and  erosion.  The 
significant  role  they  do  play  is  in  limiting  the  carburizing  potential  of 
the  truly  reactive  components  in  the  mixture.  Argon  and  Helium  have  been 
excluded  from  the  table  since  they  truly  are  "inert".  The  ratio  of  specific 
heats,  y,  has  been  included  in  the  table  at  both  ambient  and  peak  temperature 
for  each  test  gas  mixture,  to  indicate  the  changing  heat  transfer  capability 
of  a  gas  mixture  during  adiabatic  compression.  The  significance  of  convec¬ 
tive  heating  to  the  overall  results  of  this  program  will  be  discussed  in 
Section  V. 
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TABLE  3.  TEST  GAS  MIXTURE  COMPOSITION 
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V.  CORRELATION  OF  RESULTS 


5 . 1 _ Erosion  and  Flow  Conditions 


5.1.1  Pressure,  Temperature  and  Carburization  Effects 

Five  groups  of  tests  were  conducted  to  quantify  erosion  due  to 
carburization  of  4340  steel  as  described  previously  in  Section  4.2.  The  pri¬ 
mary  variable  among  groups,  the  fourth  group  being  an  exception,  was  carbon 
monoxide  concentration.  The  primary  variable  among  group  members  was  compres¬ 
sion  ratio  which  determined  pressure  and  temperature.  A  baseline  test  mixture, 
labeled  "inert",  comprised  of  45.5  percent  nitrogen  and  54.5  percent  argon 
was  selected  and  utilized  during  the  previous  two  programs. 2> 3  phis  baseline 
mixture  was  modified  in  this  program  to  include  carbon  monoxide  by  replacing 
nitrogen  with  CO.  With  the  exception  of  the  first  group  of  tests,  runs  63 
through  71,  the  total  concentration  of  carbon  monoxide  and  nitrogen  remained  ' 
constant  to  maintain  comparable  gas  temperature  and  pressure  conditions 
throughout  the  program.  In  carbon  monoxide  tests,  CO  concentration  was  kept 
relatively  high,  from  10  to  45.5  percent,  to  simulate  the  expected  CO  levels 
in  gas  evolved  from  burning  low  flame  temperature  propellants,  particularly 
nitramines.  Pressure  levels  ranged  from  193  to  335  MPa  with  the  majority 
of  tests  recording  peak  pressures  from  220  to  300  MPa. 

The  variation  of  mass  loss  with  flow  conditions,  as  specified  by 
peak  pressure  and  by  peak  temperature  are  shown  in  Figures  6  and  7,  respectively. 
The  threshold  at  which  sample  mass  loss  is  first  observed  occurs  at  248  MPa, 

3530  K  for  an  inert  (45.5%  N2)  gas  mixture,  216  MPa,  3340°K  for  a  mixture 
containing  45.5  percent  carbon  monoxide,  and  193  MPa,  3230°K  for  this  same 
45.5  percent  CO  mixture  when  4%  of  the  argon  is  replaced  with  helium.  The 
mass  loss  curve  for  this  last  mixture  is  indicated  by  two  dashed  lines  since 
it  is  unclear  from  the  point  scatter  on  this  group  of  tests,  what  the  true 
shape  of  the  curve  should  be.  (Correlation  of  heating  data  implies  that  the 
right  hand  curve  is  more  nearly  correct.)  However  the  offset  in  erosion  thresh¬ 
old  between  45.5  percent  carbon  monoxide  mixtures  that  differ  only  in  total  gas 
mass,  produced  by  substituting  helium  for  part  of  the  argon,  implies  that  the 
carbon  monoxide  surface  reaction  is  dependent  more  on  diffusion  time,  than  on 
the  surface  concentration  gradient  of  the  diffusing  species,  the  predominant 
control  in  the  case  with  oxygen.  2  Mixtures  containing  less  than  45.5  per¬ 
cent  CO,  namely  runs  83  through  88,  do  fall  between  the  45.5  percent  CO  and 
inert  gas  curves  which  would  indicate  some  diffusion  control-  at  above  threshold 
conditions,  but  since  the  measured  erosion  of  these  mixtures  does  not  approach 
the  inert  value  with  decreasing  CO  concentration  in  any  logical  fashion,  no 
true  implication  of  diffusion  controlled  erosion  can  be  made. 


Erosion  curves  for  1  and  2.5  percent  oxygen  mixtures  from  the 
previous  program2  have  been  included  on  Figures  6  and  7  for  comparison  with 
the  new  curves  representing  neutral  and  carburizing  activity.  As  was 
previously  noted2  for  the  2.5  percent  oxygen  curve,  the  shift  in  the  erosion 
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Figure  6  MASS  LOSS  vs  PEAK  TEST  GAS  PRESSURE  FOR 
VARIOUS  TEST  GAS  MIXTURES 
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Figure  7  MASS  LOSS  vs  PEAK  TEST  GAS  TEMPERATURES  FOR 
VARIOUS  TEST  GAS  MIXTURES 
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threshold  for  the  45.5  percent  carbon  monoxide  mixture  from  the  inert 
threshold  conditions  is  an  indication  that  melting  is  not  a  prerequisite 
for  the  onset  of  a  surface  reaction,  which  is  the  mechanism  by  which  carbu¬ 
rization  of  steel  enhances  material  loss.  The  parallel  nature  of  the  CO 
and  inert  curves  above  their  respective  erosion  thresholds  further  indicates 
that  the  effects  of  carburization  represent  a  nearly  constant  increment  to 
the  erosion  caused  by  convective  heating  of  4340  steel  by  propellant  gas  flow 
conditions  found  in  large  caliber  guns. 


5.1.2  Convective  and  Chemical  Heating  Effects 

Part  of  the  enhanced  CO  erosion  indicated  by  Figures  6  and  7  may 
be  attributed  to  differences  in  the  ballistic  cycle.  The  amount  of  material 
removal  is  basically  a  function  of  the  steel  sample's  entire  heating  history. 

Both  the  time  at  a  high  heating  rate  and  the  dynamic  heat  conduction  pro¬ 
cesses  are  important.  The  measured  total  heat  input  represents  the  integrated 
convective  and  chemical  heat  flux  to  the  sample  and  is  an  indicator  of 
differences  in  ballistic  cycle  time  between  tests  with  similar  heat  flux 
levels  and  chemical  environments. 

The  measured  total  heat  input,  with  the  correction  for  radial 
heat  dispersion,  as  given  in  Table  2,  is  plotted  against  mass  loss  in  Figure  8. 

A  gross  comparison  of  test  results  for  mixtures  experiencing  similar  peak 
flow  conditions  reveals  that  a  mixture  containing  CO  produces  higher  specimen 
mass  loss  at  a  lower  experimental  total  heat  input  than  a  comparable  inert 
gas  mixture.  However,  there  appears  to  be  no  quantitative  correlation  bet¬ 
ween  mass  loss  and  experimental  heat  input  in  Figure  8  because  the  act  of 
material  removal,  which  varies  in  magnitude  from  test  to  test,  removes  a 
likewise  varying  amount  of  heat  from  the  location  of  the  in-wall  thermocouple. 
Therefore,  a  correction  factor  was  derived  that  when  added  to  the  heat  input 
data  used  in  Figure  8,  would  better  correlate  mass  loss  measurements.  The 
factor  is  based  on  the  assumptions  that  any  material  removed  from  a 
sample's  surface  was  done  so  at  1720°K,  the  solidus  temperature  of  4340 
steel,  and  that  the  material  was  removed  before  its  latent  heat  was 
absorbed  by  the  remaining  steel .  The  correction  factor  is 

AQ  =  pcAxAT  =  3.27A x  J/mm2 

where  p  is  the  density  of  4340  steel  (7840  kg/m3);  c  is  the  steel's  specific 
heat  (586  J/kg-°K);  Ax  represents  the  diametral  erosion  in  millimeters  at 
the  location  of  the  heat  sensor  (column  B  of  Table  2);  and  AT  is  the 
difference  in  ambient  and  solidus  temperatures  for  the  steel  (1425°K).  A 
factor  of  1/2  is  included  in  the  equation  to  convert  diametral  recession  to 
an  averaged  radial  recession.  The  calculated  AQ  for  each  test  is  simply  added 
to  the  heat  input  given  in  Table  2  to  form  Qcorr.  Thus,  Qcorr  represents  an 
effective  heat  input  consisting  of  convective  and  chemical  components.  This 
quantity  treats  a  reduction  in  surface  temperature  required  for  melting  as  an 
equivalent  increase  in  heating.  Thus,  differences  in  material  loss  caused  by 
surface  chemistry  cause  differences  to  appear  in  Q  between  tests  with  similar 
convective  heating,  temperature,  and  pressure.  corr 
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Figure  8  MASS  LOSS  vs  TOTAL  HEAT  INPUT  FOR  TEST  GAS  MIXTURES 
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The  correlation  of  mass  loss  with  the  total  heat  input,  corrected 
for  material  loss  is  shown  in  Figure  9.  The  improved  agreement  among  tests 
in  each  group  is  readily  apparent  as  well  as  a  noticeable  erosion  path  for 
all  the  tests  for  which  heating  data  was  available.  Instrumentation  failure 
during  the  tests  for  45.5  percent  CO  mixtures  at  threshold  conditions  prevents 
the  illustration  of  a  critical  heat  input  for  this  mixture.  However,  the 
points  representing  the  45.5%  CO-4%  He-50.5%  Ar  test  group  with  the  same 
active  chemistry  as  the  45.5%  CO-54. 5%  Ar  mixture,  indicates  that  this 
carburizing  mixture  requires  less  total  heat  input,  up  to  .15  J/mm^  less  in 
some  cases,  to  affect  the  same  material  removal  than  the  "inert"  gas  mix¬ 
ture  requires.  This  variation  in  total  heating  to  produce  similar  levels 
of  erosion  can  be  attributed  to  surface  chemistry,  convective  heat  transfer 
coefficient,  heat  required  to  remove  material,  and  heat  absorption  characteris¬ 
tics.  In  order  to  delineate  the  contribution  of  each  mechanism  to  the 
heating  cycle  of  a  particular  gas  mixture  tested  in  the  STG,  a  calculation 
of  convective  hot  wall  heat  flux  for  each  test  was  made. 

Convective  hot  wall  heat  flux  is  the  driving  potential  for  melting 
erosion  once  the  melting  temperature  is  achieved.  Similarity  of  heating 
cycles  is  important  because  the  heat  flux  history  determines  the  surface 
temperature.  Therefore,  if  the  relative  difference  in  heating  cycles  is 
known,  as  is  the  case  with  45.5  percent  carbon  monoxide  mixtures  (with  and 
without  a  4  percent  substitution  of  helium  for  argon),  the  computed 
convective  heat  flux  to  the  hot  sample  surface  will  quantify  that  difference 
in  heating  cycles  since  the  heating  due  to  chemistry  remains  constant.  In 
the  opposite  case  with  test  mixtures  comprised  of  54.5  percent  argon  and 
the  remaining  mole  fraction  being  various  ratios  of  carbon  monoxide  and 
nitrogen,  convective  hot  wall  heat  flux  calculations  reflect  the  differences 
in  chemical  heating  among  tests,  namely  the  fluctuation  of  critical  surface 
temperature  at  which  simple  melting  erosion  may  commence. 

To  perform  the  heat  flux  calculation,  flow  conditions  for  each 
test  were  computed,  using  a  combination  of  two  computer  programs  and  the 
experimentally  measured  peak  pressures.  The  STG  cycle  model  (Appendix  I) 
was  used  to  create  the  proper  pressure  profile.  Also  included  in  this  par¬ 
ticular  code  is  a  representation  of  the  flow  through  the  test  sample  and 
calculation  of  the  convective  heat  flux  using  the  empirical  equation  for 
turbulent  flow  over  a  flat  plate.  The  gas  conditions  of  temperature,  density 
and  velocity  are  evaluated  in  this  code.  The  local  hot  wall  heat  flux  is 
computed  and  a  running  summation  of  the  heat  input  to  the  wall  is  evaluated. 

The  computed  total  heat  input  was  compared  with  the  experimentally  deter¬ 
mined  value  for  the  input  gas  mixture.  A  factor  was  applied  to  the  heat  flux 
calculation  to  bring  the  total  heat  input  into  agreement  with  the  experimentally 
determined  value.  In  this  manner,  the  instantaneous  value  of  convective 
heating,  as  determined  by  the  flow  conditions  and  exclusive  of  chemical 
heating,  is  believed  to  be  reasonably  correct  for  tests  with  inert  and 
carbon  monoxide  gas  mixtures. 

The  gas  temperature,  as  determined  by  the  STG  model  in  its  current 
state  of  development,  is  only  an  approximate  calculation.  A  more  accurate 
temperature  calculation  is  provided  by  the  equilibrium  combustion  code  that 
is  described  briefly  in  Appendix  II.  The  code  computes  an  isentropic 


30 


MASS  LOSS  -  mg 


800 


700 


600 


500 


400 


300 


200 


100 


0 


i 

i 

i 

! 

i 

t 

i 

t 

i 

t 

. 1 . 1 . ! . ] . 

j  !  j 

'Si! 

O  1 

0%  CO 

",  0-20°/ 

1 - 

&n2 

j 

1 

< 

i 

i 

_ 1 _ 

» 

i 

i 

i 

j 

1 

i 

i 

i 

_} 

:  i 

i  i 

i  i  i  i 

^  ..pT  ....... 

iii' 

u  lu/o  uu,  u/o  n2,  o/o  r 
□  10%  l\l2  (INERT) 
-■  45.5%  N2  (INERT) 
^  40.5-15.5%  CO,  5-30°/ 

i...  ac  co/  nr\ 

J2 

S 

.  i 

< 

i 

i 

i 

i 

_ _ i _ 

i  j  t  j 

i  i  i  ; 

i  i  i 

n2 

i 

i 

• 

i 

• 

_ « _ 

1 

— 

CD 

a  — 

■--- - 

W  .u  /O  V. 

♦  45.5%  C 

i 

_  L  _ 

:o  (4% 

He) 

|  o 

1 

68  1 

.  |  -  i  I 

i  :  i  7i 

. i . j . <_ . i  □ 

1 

1 

1 

t 

_ _ _ 

l 

1 

t 

!  •  '  1 

•  ll; 

'ii! 

1  |  i  j 

!  j  {  i 

— 

1  0 
;  66 

i 

• 

i 

i 

i 

i 

i  :  i  : 

■  i  i  i 

'  !  ' 

j  650 

1 

j 

i  i  i  i 

|  i  1  j 

i 

i 

i 

i 

l 

i 

l 

!  j  ;  } 

!  ;  i  ! 

i 

i 

> 

! 

j 

i  i  i 

i  i  i  i 

i 

i 

i 

t 

1 

i  1  1  i  - 

i!;i 

■  L  I  • 

i 

i 

* 

1  ! .  "T'" 

i  :  ;  i 

!  ;  ! 

0  67 

■  77  1 

1 

I 

1 

1 

i 

t 

. 1 . j - 1 - !• . -< 

I  M  I 

87 

00 

w« 

1  % 

1 

1 

1 

1 

:  '  :  i 

iii 

*  !  I  ' 

86  - 

90 

. . 

i 

i 

i 

i 

i 

i 

f  ‘ 

< 

i 

i 

i 

i 

_ '  __ 

1  ‘  !  ““T . 

1  1 
!  ! 
i  :  ! 

ill! 

85 

-.89 

79  ! 

i 

_  .  « _ 

1 

1 

1 

1 

1 . f . f . 1 . 

!  •  I 

1  I 

i  i  [  i 

_ 

(DA 

Mf 

I _ 

i 

■ 

75  : 

i 

i 

t 

i-  1  !  1 

iii! 

i 

i 

i 

i 

i 

i 

f 

. r~~i . f . - . 

i  i 

iii 

iii 

k. _ 

94 

-w-4 

74  7: 

t 

1 

l 

MZ?., 

i 

i 

i 

« 

i 

i 

{ 

. f . - . j . t . 

iii 

:  :  i  i_ 

93 

i 

» 

_ ! _ 

1 

2 

_ a 

i 

i 

i 

i 

i 

i 

« 

i 

i 

i 

i 

r — y . . i . 

ii|! 

- : - ! - 1  I _ 

1-00  1.10  1.20  1.30  1.40  1.50  1.60  1.70 

CORRECTED  TOTAL  HEAT  INPUT  (QC0RR)  -  J/mm2 

Figure  9  MASS  LOSS  vs  TOTAL  HEAT  INPUT  CORRECTED  FOR  LOCAL  HEAT 
LOST  BY  MATERIAL  REMOVAL 


31 


compression  of  gases  beginning  with  an  arbitrary  mixture.  In  addition  to 
determining  the  temperature  and  pressure,  the  concentration  of  the  various 
chemical  species  formed  during  the  equilibrium  combustion  process  are  also 
determined.  For  example,  the  mole  fraction  of  cyanide  formed  by  compression 
of  gas  mixtures  containing  nitrogen  and  carbon  monoxide  is  evaluated. 


An  approximate  technique  was  used  to  correct  the  heat  flux  calculated 
by  the  STG  model  to  the  more  accurate  temperature  conditions  of  the  equili¬ 
brium  combustion  code.  The  heat  flux  to  a  surface  is  equal  to  the  product  of 
a  coefficient  and  the  temperature  difference  between  the  gas  and  the  surface, 

q  =  h  (To  "  Tw)  • 


For  the  heat  flux  to  a  flat  plate  in  turbulent  flow,  the  coefficient  is 
functionally  proportional  to  the  gas  density,  velocity  and  viscosity, 

h  ~  (pu)0,8  y0,2 

These  gas  parameters  can  be  expressed  in  terms  of  temperature  as  follows: 

p  ~  T"1  through  the  ideal  gas  equation  of  state 
u  ~  T^<8  through  the  energy  equation, 
and  y  ~  TO • 5  from  molecular  transport  theory. 

Thus,  the  approximate  dependence  of  h  on  the  gas  temperature  is 
h  ~  (T“l  •  (pO .  5-j  0 . 2  =  j- .  3 

The  convective  hot  wall  heat  flux  computed  by  the  STG  code  is 

qS  ~  TgS  (Tgs  "  Tws) 

where  Tc,s  and  Tws  are  the  respective  peak  values  for  gas  and  surface  tempera¬ 
tures.  Similarly,  the  convective  heat  flux  to  a  melting  steel  surface  is 

%w  ~  Tge  ®  ‘  ^  (Tge  -  Twm) 


where  Tge  is  the  gas  temperature  calculated  by  the  equilibrium  combustion  code 
and  Twm  is  the  solidus  temperature  of  4340  steel,  1720°K.  The  corrected  value 
for  the  heat  flux,  qhw>  is 


This  represents  an  approximate  value  of  convective  heat  flux  to  a  melting 
surface . 
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A  graph  showing  the  correlation  of  mass  loss  with  convective  hot 
wall  heating  is  shown  in  Figure  10.  The  expected  offset  between  flux  curves 
for  45.5  percent  CO  mixtures  both  with  and  without  a  4  percent  substitution 
of  helium  for  argon  represents  the  difference  in  heat  cycle  histories  for 
gas  mixtures  of  equal  chemical  activity,  but  different  total  mass. 

The  test  gas  mixtures  of  45.5  percent  carbon  monoxide  and  45.5 
percent  nitrogen,  both  containing  identical  mole  fractions  of  argon,  have 
nearly  identical  masses,  ratios  of  specific  heat  and,  therefore,  similar 
heating  cycles.  However  the  CO  mixture  was  capable  of  producing  similar 
erosion  to  the  "inert"  gas  mixture  at  a  heat  flux  rate  nearly  100  J/mm^-sec 
lower  than  the  "inert"  mix,  both  at  and  well  above  erosion  threshold  condi¬ 
tions.  To  account  for  the  difference  in  erosion  threshold,  the  above  convec¬ 
tive  heat  flux  equation  was  recalculated  for  a  reduced  hot  wall  melting 
temperature  that  would  shift  the  CO  erosion  curve  to  coincide  with  the  "inert" 
erosion  curve.  The  computed  solidus  temperature,  assuming  all  other  proper¬ 
ties  of  the  steel  remained  the  same,  was  lowered  to  957° K.  A  reduction  of 
over  750°K  in  melting  temperature  is  unrealistic  for  simple  diffusion  of 
carbon  into  steel,  according  to  the  iron-carbon  phase  diagram  (Figure  28). 
Therefore  the  shift  in  the  CO  erosion  threshold  suggests  the  strong  influence 
of  surface  phenomena  that  alter  the  thermal  properties  of  the  materials, 
enhance  the  net  heat  input,  or  both. 


5.2  Erosion  Distribution 


Measurements  of  diametral  recession  were  made  at  the  four  axial 
locations  identified  in  Figure  3  and  are  presented  for  each  test  sample  in 
Table  2.  While  examining  these  measurements,  it  was  observed  that  the  dis¬ 
tribution  of  material  erosion  was  altered  by  the  gas  mixture.  Several  of 
these  distributions,  normalized  to  the  most  forward  measurement,  are  shown 
plotted  in  Figure  11.  In  this  figure,  the  averages  for  several  runs  in  a 
group  with  similar  characteristics  are  compared. 

The  typical  inert  distribution,  as  illustrated  by  averages  of 
points  from  Run  75  and  77,  shows  the  highest  erosion  near  the  leading  edge 
of  the  sample  which  progressively  decreases  until  very  little  or  no  erosion 
is  observed  by  the  third  measuring  point,  which  is  labeled  point  C.  Thus, 
with  only  argon  and  nitrogen  in  the  gas  mixture,  negligible  erosion  is 
observed  over  the  latter  half  of  the  sample. 

The  downstream  level  of  erosion  tends  to  increase  when  carbon 
monoxide  is  added  to  the  gas  mixture.  Run  Numbers  78  through  80  and  83  and 
84  all  have  quite  similar  gas  temperatures  and  pressures  and  all  experience 
measurable  mass  loss  during  the  test.  The  significance  of  the  erosion  dis¬ 
tribution  from  these  tests  in  comparison  with  that  of  the  inert  test,  is 
that  the  erosion  at  Points  B  and  C  is  substantially  higher. 
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Figure  10  INFLUENCE  OF  CONVECTIVE  HEATING  AND  CARBURIZING-TYPE 
REACTIONS  ON  EROSION  OF  4340  STEEL 
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The  average  of  Runs  85  through  88  show  an  overwhelming  increase 
in  erosion  towards  the  downstream  end  of  the  test  sample.  The  erosion  of 
the  downstream  end  of  the  sample  in  relation  to  that  of  the  leading  edge, 
which  was  used  to  normalize  results,  is  likely  not  as  dramatic  as  indicated 
in  Figure  11  because  the  erosion  at  the  leading  edge  appears  to  be  unusually 
small.  However,  erosion  at  the  downstream  end  of  the  samples  is  large  in 
relation  to  that  of  the  inert  samples  and,  therefore,  it  appears  as  though 
carbon  monoxide  in  the  gas  is  contributing  to  a  downstream  erosion  effect. 

This  is  similar  to  that  observed  when  oxygen  was  introduced  into  the  test 
gas  as  described  in  Reference  2. 

Carbon  monoxide  and  nitrogen  have  the  same  molecular  weight.  When 
one  is  exchanged  for  the  other  in  the  test  gas,  there  is  very  little  effect 
on  STG  test  conditions.  Carbon  monoxide  can  be  substituted  for  the  nitrogen 
and  the  effects  of  carbon  monoxide  on  the  specimen  can  be  observed  with 
little  change  in  convective  heating  and  cycle  time  as  long  as  the  concentra¬ 
tion  of  argon  remains  the  same.  Therefore,  differences  in  these  quantities 
are  not  responsible  for  the  observed  differences  in  erosion. 

Rapid  falloff  of  erosion  was  observed  with  no  carbon  monoxide  in 
the  test  gas.  When  45.5%  carbon  monoxide  and  no  nitrogen  was  incorporated, 
the  falloff  was  less  rapid  as  shown  in  Figure  11,  but  very  little  or  no 
erosion  was  observed  again  at  the  trailing  edge.  When  carbon  monoxide  and 
nitrogen  were  approximately  equal  in  concentration,  that  is,  when  15  to 
30%  of  the  gas  was  carbon  monoxide,  the  maximum  enhancement  of  erosion  distri¬ 
bution  was  observed.  This  suggests  that  the  presence  of  nitrogen  in  the 
gas  plays  a  role  in  the  observed  erosion  pattern.  Again  it  is  noted  that  the 
low  erosion  at  the  leading  edge  of  the  sample  during  these  tests  caused  the 
curve  to  exaggerate  this  effect,  but  enhanced  downstream  erosion  does  appear 
to  be  real . 

This  is  also  substantiated  by  the  results  of  Runs  63  and  71  which 
are  not  shown  in  Figure  11.  Run  63  was  a  test  with  10%  carbon  monoxide  and 
90%  argon.  This  test  was  very  high  in  temperature  and  a  great  deal  of  material 
was  lost.  The  diametral  recession  was  nearly  constant  along  the  entire 
length  of  the  sample  with  the  measurements  being  equal  at  both  the  leading 
and  trailing  edges.  Run  71  was  a  test  with  10%  nitrogen,  no  carbon  monoxide, 
and  90%  argon.  The  pressure  and  temperature  of  these  two  runs  was  almost 
identical.  The  erosion  at  the  leading  edge  of  Run  71  was  the  same  as  for 
Run  63.  However,  at  the  trailing  edge,  erosion  was  approximately  half  that 
of  the  forward  measurement. 

It  is  evident  that  surface  chemistry  plays  an  important  role  in  mate¬ 
rial  loss  in  a  carbon  monoxide  atmosphere.  The  inert  run  indicates  a  falling 
off  of  erosion  in  much  the  way  that  one  would  expect  heat  transfer  to  fall 
off  in  turbulent  flow  over  a  flat  plate.  When  carbon  monoxide  is  present, 
the  falloff  in  erosion  is  altered.  The  data  show  the  observed  effect  to 
be  independent  of  convective  heating.  The  resemblance  of  this  effect  to 
that  of  oxygen  cited  earlier  suggests  the  presence  of  an  exothermic  reaction. 
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NORMALIZED  DIAMETRAL  RECESSION 


Figure  11  INFLUENCE  OF  TEST  GAS  CONSTITUENTS  ON  SURFACE  RECESSION 
OF  4340  STEEL 
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5.3 


Surface  Characterization 


The  Scanning  Electron  Microscope  was  used  to  help  characterize  the 
surface  features  of  each  test  sample.  Selected  photographs  are  shown  in 
subsequent  figures.  The  first  four  figures  (12-15)  illustrate  the  surface 
of  an  inert  sample  with  increasingly  severe  test  conditions.  Figures  12 
and  13  show  the  sample’s  surface  at  the  inlet,  center,  and  exit  point  for 
flow  conditions  at  the  onset  of  melting.  That  is,  the  point  where  weight 
loss  is  first  observed.  At  this  time,  erosion  is  initiated  at  the  leading 
edge  of  the  sample.  As  conditions  increase,  the  point  at  which  erosion 
occurs  moves  farther  downstream  on  the  sample.  These  figures  show  the 
presence  of  small  solidified  droplets.  There  is  also  what  appears  to  be 
a  small  crust  of  solidified  material..  The  presence  of  a  crust- like  layer 
is  more  evident  in  Figures  14  and  15  which  are  representative  of  samples 
tested  in  an  inert  atmosphere  above  the  erosion  threshold.  The  surface  of 
these  samples  displays  a  mottled  characteristic  that  appears  to  be  solidified 
melt  or  surface  layer  of  some  sort. 

The  next  series  of  four  figures  (16-19)  illustrates  surface  features 
of  4340  steel  samples  tested  in  an  atmosphere  consisting  of  45.5%  carbon 
monoxide  and  argon.  The  difference  between  these  and  the  inert  tests  was 
that  the  nitrogen  was  exchanged  for  carbon  monoxide.  The  molecular  weight  of 
these  two  gases  is  the  same  and  test  conditions  and,  most  significantly,  the 
convective  heating  are  also  approximately  the  same.  The  samples  shown  in 
Figures  16  and  17  were  tested  at  approximately  the  same  conditions  and  show 
what  appear  to  be  solidified  droplets  or  perhaps  molten  protrusions  that 
existed  on  the  base  metal  surface.  Otherwise,  there  is  very  little  change 
to  the  initial  surface.  At  more  severe  conditions  (Figures  18  and  19)  it  is 
noticed  that  the  surface  appears  to  be  quite  clean  although  rough  and  exhibits 
what  might  be  interpreted  as  a  turbulent  flow  pattern.  The  surface  does  not 
have  the  crust-like  layer  that  was  observed  for  the  inert  shots.  One  might 
postulate  that  the  presence  of  carbon  monoxide  facilitated  removal  of  the 
crust- like  layer  that  was  observed  in  the  inert  shots  and  that  this  is 
also  responsible  for  the  increased  erosion  at  farther  downstream  stations 
that  was  noted  in  the  previous  section. 

Some  cracking  was  also  noted  in  the  regions  of  high  erosion.  This 
is  contrary  to  the  observation  made  in  Reference  2,  concerning  the  effects 
of  oxidation,  where  significant  cracking  was  observed  at  the  onset  of  erosion 
rather  than  in  a  region  of  high  erosion. 

The  next  series  of  four  figures  (20-23)  illustrates  the  surface 
features  of  the  4340  carbon  steel  sample  as  influenced  by  decreasing  percen¬ 
tage  of  carbon  monoxide.  For  these  runs,  the  test  conditions  are  very  much 
the  same  and  nitrogen  was  exchanged  for  carbon  monoxide  in  an  attempt  to  fill 
in  the  region  between  inert  and  45.5%  carbon  monoxide.  The  significant  feature 
of  these  SEM  photographs  is  the  slight  tendency  of  the  surface  features 
to  approach  those  of  an  inert  sample  with  decreasing  carbon  monoxide  concen¬ 
tration.  When  the  gas  contains  15.5%  carbon  monoxide,  the  lowest  concentration 
tested,  the  surface  still  tends  to  resemble  more  of  the  high  concentration  of 
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carbon  monoxide  rather  than  the  inert  sample.  This  seems  to  indicate  that 
the  process  is  only  slightly  diffusion  limited  for  the  concentration  tested. 
An  interesting  aspect  of  the  40.5%  carbon  monoxide  run  is  that  it  was  tested 
at  a  higher  pressure  than  any  of  the  others  and  the  flow  patterns  that  were 
observed  looked  like  eyes  on  a  potato.  This  indicates  a  high  interaction  of 
the  turbulent  flow  with  the  surface  and  the  presence  of  a  sluffing  mechanism. 
At  the  lower  concentrations  of  carbon  monoxide,  this  feature  was  not  apparent 
At  low  to  intermediate  concentrations  of  carbon  monoxide,  the  presence  of 
solidified  melt  became  more  evident,  whereas  at  the  high  concentrations  all 
the  melt  appeared  to  be  swept  away. 

The  final  set  of  four  figures  (24-27)  illustrates  surface  features 
of  a  sample  tested  in  a  gas  mixture  containing  4%  helium,  where  helium  was 
exchanged  for  some  of  the  argon.  This  series  of  four  tests  had  gradually 
increasing  pressures.  Substantial  erosion  was  observed  at  the  three  higher 
pressures  where  201  MPa  appears  to  be  the  threshold  point.  The  surface 
features  of  all  three  appear  to  resemble  each  other  with  respect  to  the 
formation  of  a  kind  of  a  web.  A  dendrite  feature  on  the  surface  could  be 
solidification  of  a  surface  melt  or  a  slushy  layer.  Flow  turbulence  near 
the  surface  and  some  carcking  was  observed  under  all  conditions.  The  marbled 
crust  that  appeared  with  the  inert  shot  was  not  present  here. 


5.4  Summary  of  Results 

Carbon  monoxide  gas  interacts  chemically  with  a  gun  barrel  wall 
within  the  time  constraints  of  a  ballistic  cycle.  Evidence  of  surface  phe¬ 
nomena  between  CO  and  steel  to  support  this  statement  are: 

1*  A  shift  in  the  onset  of  erosion  due  to  the  presence 
of  carbon  monoxide. 

2.  Enhancement  of  downstream  erosion  when  carbon  mono¬ 
xide  is  a  test  gas  constituent. 

3.  A  change  in  surface  characteristics  of  barrel 
steel  when  exposed  to  carbon  monoxide  rather 
than  to  nitrogen. 

Several  explanations  for  the  short  time  CO-steel  erosion  mechanism 
are  postulated.  Simple  carburization  or  diffusion  of  carbon  atoms,  into 
the  steel  surface  would  produce  several  effects.  First,  increased  carbon 
concentration,  present  either  as  graphite  or  iron  carbide,  would  reduce 
the  solidus  temperature  of  the  steel  as  shown  on  the  Iron-Carbon  Equilibrium 
Diagram  in  Figure  28.  Second,  an  increase  in  surface  carbon  concentration 
would  reduce  the  thermal  conductivity  of  the  affected  layer,  causing  it  to 
experience  an  accelerated  temperature  rise.  The  elevated  heat  content  of 
the  affected  layer  would  be  further  enhanced  by  the  exothermic  iron  carbide 
(cementite)  reaction  between  CO  and  iron,  even  when  flow  conditions  started 
to  become  less  severe. 
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A  second  explanation  of  the  CO-steel  erosion  mechanism,  previously 
proposed, 4  attributes  increased  heat  input  to  the  entire  steel  surface  layer 
from  the  highly  exothermic  boundary  layer  reaction  of  carbon  monoxide  molecules 
and  iron  to  form  iron  penta  carbonyl,  FeCCO)^.  To  a  lesser  degree,  carbon 
monoxide  molecules  may  also  be  reacting  with  chromium  and  nickel,  the  chief 
alloying  metals  in  barrel  steel  to  produce  equally  volatile  Cr(CO)^  and 
Ni (C0)4 . 


At  the  present  time,  metal lographic  examination  of  the  4340  steel 
samples,  tested  in  this  program,  are  being  conducted  at  Government  laboratories. 
It  is  anticipated  that  these  additional  examinations  will  shed  light  on  the 
precise  erosion  mechanisms  that  can  be  linked  to  propellant  gas  atmospheres, 
rich  in  carbon  monoxide. 
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SAMPLE  38,  RUN  72,  1000  MAGNIFICATION 

4340  STEEL  SAMPLE  SURFACE  TESTED  IN  AN 
"INERT"  ATMOSPHERE  TO  223  MPa,  3422°K 
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FLOW 
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EXIT 


SAMPLE  39,  RUN  73,  1000  MAGNIFICATION 

Figure  13  4340  STEEL  SAMPLE  SURFACE  TESTED  IN  AN 
"INERT"  ATMOSPHERE  TO  248  MPa,  3531°K 
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EXIT 


SAMPLE  41,  RUN  75,  1000  MAGNIFICATION 

Figure  14  4340  STEEL  SAMPLE  SURFACE  TESTED  IN  AN 
"INERT"  ATMOSPHERE  TO  285  MPa,  3678°K 


FLOW 
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Figure  15  4340  STEEL  SAMPLE  SURFACE  TESTED  IN  AN 
"INERT"  ATMOSPHERE  TO  335  MPa,  3865°K 
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EXIT 


SAMPLE  45,  RUN  81,  1000  MAGNIFICATION 


Figure  16  4340  STEEL  SAMPLE  SURFACE  TESTED  IN  A  45.5  PERCENT 
CARBON  MONOXIDE  ATMOSPHERE  TO  213  MPa,  3322°K 
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INLET 


SAMPLE  40A,  RUN  82,  1000  MAGNIFICATION 

Figure  17  4340  STEEL  SAMPLE  SURFACE  TESTED  IN  A  45.5  PERCENT 
CARBON  MONOXIDE  ATMOSPHERE  TO  216  MPa,  3338°K 
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SAMPLE  43,  RUNS  78,  79,  1000  MAGNIFICATION 

Figure  18  4340  STEEL  SAMPLE  SURFACE  TESTED  IN  A  45.5  PERCENT 
CARBON  MONOXIDE  ATMOSPHERE  TO  277  MPa,  3588°K 
AND  256  MPa,  3507°K 
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EXIT 


SAMPLE  44,  RUN  80,  1000  MAGNIFICATION 


Figure  19  4340  STEEL  SAMPLE  SURFACE  TESTED  IN  A  45.5  PERCENT 
CARBON  MONOXIDE  ATMOSPHERE  TO  245  MPa,  3462°K 
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Figure  20  4340  STEEL  SAMPLE  SURFACE  TESTED  IN  A  40.5  PERCENT 
CARBON  MONOXIDE  ATMOSPHERE  TO  299  MPa,  3690°K 
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SAMPLE  49,  RUN  86, 1000  MAGNIFICATION 

Figure  21  4340  STEEL  SAMPLE  SURFACE  TESTED  IN  A  25.5  PERCENT 

CARBON  MONOXIDE  ATMOSPHERE  TO  281  MPa,  3637°K 
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EXIT 


INLET 


SAMPLE  50,  RUN  87,  1000  MAGNIFICATION 


Figure  22  4340  STEEL  SAMPLE  SURFACE  TESTED  IN  A  20  5  PERCENT 
CARBON  MONOXIDE  ATMOSPHERE  TO  287  MPa,  36650K 
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EXIT 


INLET 


SAMPLE  51,  RUN  88,  1000  MAGNIFICATION 

Figure  23  4340  STEEL  SAMPLE  SURFACE  TESTED  IN  A  15.5  PERCENT 
CARBON  MONOXIDE  ATMOSPHERE  TO  288  MPa,  3673°K 
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EXIT 


SAMPLE  52,  RUN  89,  1000  MAGNIFICATION 

Figure  24  4340  STEEL  SAMPLE  SURFACE  TESTED  IN  A  45.5  PERCENT 
CARBON  MONOXIDE  ATMOSPHERE  TO  250  MPa,  3479°K 
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INLET 


SAMPLE  54,  RUN  91,  1000  MAGNIFICATION 

Figure  25  4340  STEEL  SAMPLE  SURFACE  TESTED  IN  A  45.5  PERCENT 
CARBON  MONOXIDE  ATMOSPHERE  TO  236  MPa,  3422°K 
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.  SAMPLE  55,  RUN  92, 1000  MAGNIFICATION 

Figure  26  4340  STEEL  SAMPLE  SURFACE  TESTED  IN  A  45.5  PERCENT 
CARBON  MONOXIDE  ATMOSPHERE  TO  219  MPa,  3351°K 
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SAMPLE  57,  RUN  94,  1000  MAGNIFICATION 

Figure  27  4340  STEEL  SAMPLE  SURFACE  TESTED  IN  A  45.5  PERCENT 
CARBON  MONOXIDE  ATMOSPHERE  TO  201  MPa,  3268°K 
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Figure  28  IRON-CARBON  EQUILIBRIUM  DIAGRAM 
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Temperature, 


VI.  CONCLUSIONS 


This  program  is  part  of  an  ongoing  effort  to  isolate  the  various 
phenomena  that  contribute  to  gun  barrel  erosion.  Research  dealing  with 
erosion  of  4340  steel  in  a  carburizing  atmosphere,  the  primary  topic  of 
investigation,  led  to  the  following  conclusions: 

1.  Addition  of  carbon  monoxide  to  the  test  gas  mixture 
was  shown  to  increase  the  erosion  of  4340  steel. 

2.  The  increase  in  erosion  was  nearly  independent  of 
carbon  monoxide  content  which  is  indicative  of 
reaction  rate  controlled  carburization. 

3.  The  addition  of  carbon  monoxide  shifted  the  erosion 
threshold  to  less  severe  flow  conditions. 

4.  An  inert  gas  mixture  (45.5%  N2)  produced  no  indica¬ 
tion  of  surface  melting  at  conditions  comparable  to 
the  erosion  threshold  for  a  gas  mixture  containing 
45.5  percent  carbon  monoxide.  This  indicates  that 
carburization  takes  the  form  of  a  surface  reaction. 
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APPENDIX  I 


SHOCK  TUBE  GUN  COMPUTER  SIMULATION 


Overview 


The  mathematical  model  described  in  this  Appendix  simulates  the 
operation  of  the  Shock  Tube  Gun.  The  simulation  provides  a  complete  descrip¬ 
tion  of  the  entire  cycle  of  the  Shock  Tube  Gun,  beginning  with  release  and 
subsequent  acceleration  of  the  piston  by  the  high  pressure  driver  gas.  As 
the  piston  is  accelerated  through  the  driven  tube,  the  simulation  computes 
the  increase  in  pressure  and  temperature  of  the  test  gas.  In  addition,  the 
simulation  evaluates  the  total  temperature  pressure  and  density  in  the 
plenum  chamber  and  computes  the  flow  through  the  test  specimen.  The  heat 
flux  to  the  specimen,  the  resulting  temperature  history  at  a  location  on 
the  surface  of  the  specimen  and  the  temperature  distribution  normal  to 
the  surface  are  also  calculated.  Finally,  the  simulation  calculates  the 
travel  of  a  projectile  through  the  barrel. 

The  objective  of  this  code  is  to  provide  a  means  for  calculating 
test  conditions  for  the  purpose  of  establishing  the  initial  driver  pressure 
and  gas  mixture.  Differences  in  the  ballistic  cycle  of  the  Shock  Tube  Gun 
due  to  gas  composition  are  reflected  through  differences  in  pressure  and 
total  heat  input.  A  primary  use  of  the  code  is  to  help  distinguish  between 
erosion  due  to  melting  and  that  due  to  chemical  effects.  This  is  done 
through  computation  of  the  convective  heating  to  the  sample  without  chemistry 
which  provides  a  means  for  comparing  tests  within  a  test  matrix  on  an  equal 
basis  as  far  as  the  inherent  flow  heating  of  the  test  gas.  Thus,  excess 
material  removal  from  one  gas  mixture  in  comparison  to  another  is  likely 
due  to  chemical  effects.  This  further  allows  estimates  to  be  made  of  the 
effective  heat  input  due  to  chemical  effects.  This  can  be  done  by  comparing 
the  heat  input  at  the  onset  of  erosion,  or  at  points  of  equal  erosion  between 
inert  and  chemically  active  gases.  The  code  enables  the  facility  user  to 
quantify  levels  of  heating  experienced. 

The  major  assumption  applied  with  formulating  the  code  was  that 
of  quasi-steady  operation.  That  is,  pressure  waves  and  other  unsteady  aspects 
of  the  event  are  not  calculated.  The  pressure  is  assumed  to  be  constant 
throughout  the  driver  system  and  throughout  the  driven  tube  at  any  instant 
of  time  during  the  compression  cycle. 

The  other  limiting  assumption  that  is  currently  employed  In  this 
code  is  that  of  a  frozen  gas  composition  whereby  the  initial  gas  composition 
is  assumed  to  be  maintained  throughout  the  ballistic  cycle.  This  assumption 
influences  the  resulting  temperature  and  pressures  to  some  extent  in  cases  wher 
chemical  reactions  and  dissociation  become  important. 


The  individual  gas  constituents  are  assumed  to  be  mixed  uniformly. 
Furthermore,  the  parameters  used  by  the  van  der  Waals  equation  of  state  for 
an  imperfect  gas  and  the  temperature  dependency  of  specific  heat  are  assumed 
to  be  satisfied  through  a  linear  averaging  according  to  mole  fraction. 

The  assumption  of  frozen  gas  constituency  is  an  interim  assumption 
that  will  be  relieved  when  a  chemical  equilibrium  gas  code,  also  discussed 
in  this  appendix,  can  be  combined  with  the  STG  model.  At  that  time,  the  gas 
constituency  will  be  assumed  to  be  in  thermal  and  chemical  equilibrium  at 
all  times.  At  present,  the  discrepancy  between  STG  code  calculations  and 
actual  STG  test  results,  as  illustrated  in  the  figure  below,  can  be  attributed 
to  both  the  quasi-steady  limitations  of  the  code  and  the  inherent  accuracy 
limitations  of  the  test  instrumentation. 


Figure  PEAK  TEST  GAS  PRESSURE  vs  DRIVER  GAS  PRESSURE 
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Model  Description 


Piston  Motion 

The  piston  motion  is  evaluated  by  applying  a  force  balance  on  the 
piston.  The  accelerating  force  is  applied  by  the  high-pressure  nitrogen 
driver,  500  to  700  psi,  on  the  upstream  side  of  the  piston.  The  driven  gas 
on  the  downstream  side  of  the  piston  is  initially  at  atmospheric  pressure. 

It  is  assumed  that  no  gas  leaks  past  the  piston.  This  is  essentially  verified, 
at  least  initially,  by  the  maintenance  of  a  perfect  seal  and  the  ability  to 
evacuate  the  driven  gas  chamber.  The  nitrogen  gas  is  represented  by  the 
ideal  gas  equation  of  state.  The  driver  conditions  are  calculated  from 
conservation  of  energy  principles  which  are  used  to  continually  evaluate  the 
amount  of  energy  that  is  being  transferred  from  the  gas  to  the  piston.  As 
cited  previously,  this  calculation  is  quasi-steady  in  that  the  unsteady 
expansion  aspects  are  not  considered  and  the  pressure  is  assumed  to  be 
constant  throughout  the  driver  system. 

The  driver  gas  (nitrogen)  properties  are  assumed  to  be  constant 
over  the  range  of  temperature  and  pressure  encountered  during  the  compression 
cycle  and  are  specified  by: 

Equation  of  state  gas  constant,  R  =  55. 0  ft-lbf/lbm°R 

Specific  heat  at  constant  volume,  cv  =  0.177  Btu/lbm-°R 

Specific  heat  at  constant  pressure,  cp  =  0.248,  and 

Ratio  of  specific  heats,  Y  =  1.4 

The  assumption  of  a  perfect  seal  at  the  piston  infers  the  existence 
of  constant  gas  mass  in  the  driver  system  during  the  cycle  so  that 

_  V1 

PD  "  V 


where  R  and  m  are  constants,  the  driver  volume,  V,  is  expressed  in  terms  of 
piston  travel  and  the  initial  volume  by 

v  *  y  *  vy 


where  Aj  is  the  driver  tube  area. 
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Gas  temperature  is  expressed  in  terms  of  internal  energy  where 


E  =  cVmDT 

and  the  energy  change  resulting  from  work  expended  through  piston  motion 
is : 


&E  =  PDVV778'  ■  VdAT 

The  driver  gas  is  initially  at  room  temperature  and  it  is  assumed  that  the 
small  temperature  decrease  during  the  cycle  is  not  influenced  by  heat  transfer. 

Piston  motion  is  evaluated  by  applying  a  force  balance  across  the 
piston,  taking  into  account  the  frictional  drag. 


F  =  Vpd  -  Pt5  -  d 

where  pp  and  pp  are  the  respective  pressures  of  the  driver  and  test  gases,  and 
D  is  the  frictional  drag  of  the  piston  which  is  expressed  in  terms  of  piston 
velocity  by: 


Piston  acceleration,  velocity  and  travel  follow: 


F 

a  =  — 
p  m 


AV  =  anAt 
P  P 


Axp  =  -P-At2  +  v0At 


Test  Gas  Compression 

The  compression  of  the  test  gas  occurs  as  a  result  of  piston 
motion.  Energy  that  is  added  by  virtue  of  the  compression  is  calculated 
from  the  conservation  of  energy  equation.  The  work  done  by  the  piston  on  the 
test  gas  during  this  compression  is  one  term  in  this  conservation  of  energy 
equation.  Other  terms  include  heat  loss  to  the  wall  of  the  tube,  which  becomes 
important  as  the  gas  temperature  rises.  The  other  important  equation  is  con¬ 
servation  of  mass.  The  test  chamber  is  not  a  closed  chamber  but  contains 
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an  exhaust  port  at  the  downstream  end  where  the  test  sample  is  located.  Thus, 
test  gas  is  allowed  to  flow  from  the  test  chamber  through  the  test  sample. 
Therefore,  the  mass  in  the  system  is  not  constant  by  virtue  of  mass  and 
energy  flow  from  the  driven  tube  and  plenum  chamber  through  the  test  sample. 
Terms  in  the  conservation  of  energy  and  mass  equations  reflect  this  mass 
and  energy  loss. 

The  equation  of  state  that  applies  to  the  test  gas  is  the  van  der 
Waals  equation  which  includes  terms  to  express  the  nonlinear  relationship 
between  pressure  density  and  temperature.  The  terms  for  this  equation  are 
determined,  as  mentioned  previously,  by  a  linear  averaging  of  the  mole 
fraction  of  the  test  gas  constituents. 

The  test  gas  specific  heat  is  also  assumed  to  be  for  non-perfect 
gas  and  is  expressed  in  terms  of  a  linear  function  of  the  gas  temperature. 

The  coefficients  in  this  expression  are  also  linear  averages  of  the  mole 
fraction  of  the  test  gas  compositions.  The  test  gas  specific  heat  is  expressed 
in  terms  of  a  secant  function  in  which  the  product  of  the  specific  heat  and 
the  temperature  yield  the  internal  energy.  This  is  contrary  to  a  normal 
expression  of  specific  heat  whereby  it  is  a  tangent  function  so  that  the 
integral  of  the  product  of  specific  heat  and  temperature  yield  the  internal 
energy.  The  technique  used  here  provides  rather  simple  yet  effective  means 
for  evaluating  the  internal  energy  in  a  finite  difference  scheme  with 
many  time  steps  . 

The  van  der  Waals  equation  of  state  used  for  the  test  gas  is 


a. 

where  -g— 
S1 


27  ai 

RT~-  —  -  27Pci- 


8 


ci’  S  2 
1 


v  is  the  specific  volume,  and  Tci  and  pci  are  the  critical  temperature  and  pres¬ 
sure  for  the  ith  gas  constituent,  a  and  3  are  the  average  quantities  based  on 
the  mole  fraction  of  the  ith  constituent. 


The  specific  heat  at  constant  volume  is  defined  by 


c 

v 


C(c 


VO 


Z.X  . 
i  gi 


cvT.(T  -  460)ei)T 


ap  ,  1 

778J  Ti 


where  cvo  and  cvTi  are  the  intercept  and  slope  of  the  temperature-dependent 
specific  heat  and  XCTi  is  the  mole  fraction  of  the  ith  constituent. 

to 


which  includes  temperature  and  high  density  effects. 

The  test  chamber  mass  balance  is  given  by 
"V  =  mTo  -  Am 
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where  Am  is  the  mass  flow  through  the  test  sample.  Initially  this  flow  is 
assumed  to  be  negligible  and  the  test  chamber  and  barrel  volumes  are  lumped 
together.  When  the  projectile  velocity  exceeds  100  ft/sec,  the  calculation 
of  flow  through  the  test  sample  is  initiated.  This  computation  involves 
determination  whether  sonic  or  subsonic  flow  conditions  exist  within  the 
test  sample.  The  sonic  static  pressure  is  given  by: 


P* 


=  PXC1  + 


Y  -  U-Y/V  -  1 
2  j 


where  p^  is  the  test  chamber  pressure. 

If  p*  is  greater  than  the  barrel  pressure,  P2,  (downstream  from  the  test 
sample)  then  sonic  conditions  exist  and 

l /2 

Am  =  P*A12(rj5')  At, 

where  A]_2  is  the  test  sample  flow  area,  T*  is  the  sonic  static  temperature, 
and  At  is  the  computation  time  interval . 

If  subsonic  conditions  exist, 


Am  =  8.02  A12(Cpi  -  P2)  ^p)1/2At, 

where  is  the  equation  for  flow  through  a  venturi  in  terms  of  the  upstream 
and  downstream  line  pressures. 

The  change  in  internal  energy  in  the  chamber  over  the  calculation 
time  interval  is  given  by 

plVXp 

^E1  =  778  "  cvAmTlY 


which  represents  the  compression  work  due  to  the  piston  and  the  loss  in 
enthalpy  due  to  flow  through  the  test  sample.  Heat  transfer  losses  in  the 
driven  tube  and  chamber  are  not  included  in  this  analysis  at  present.  The 
gas  chamber  temperature  is  defined  by 


where  is  the  current  value  of  internal  energy. 
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The  Test  Sample 


The  test  sample  is  a  straight  channel  with  a  radiused  entrance. 

The  flow  through  the  channel  is  computed  by  either  sonic  or  subconic  condi¬ 
tions  depending  on  the  pressure  at  the  inlet  and  outlet  to  the  sample.  The 
calculations  provide  the  static  flow  conditions  of  pressure,  temperature,  and 
density  in  addition  to  the  flow  velocity  over  the  surface  of  the  sample.  These 
conditions  are  in  turn  used  in  the  equation  that  expresses  turbulent  heat 
flux  to  a  flat  plate.  The  heat  flux  is  computed  and  then  summed  to  yield  a 
current  level  of  total  heat  input.  The  heat  flux  to  the  surface  is  also 
applied  to  an  unsteady  heat  conduction  routine  by  which  the  surface  tempera¬ 
ture  and  the  temperature  distribution  in  the  test  sample  are  evaluated. 

These  calculations  are  all  performed  within  the  same  time  step  of  the  overall 
finite  difference  calculation. 

The  technique  used  to  calculate  heat  flux  to  the  test  sample 
surface  requires  the  flow  velocity,  density,  and  viscosity.  The  density  and 
viscosity  are  evaluated  on  the  basis  of  a  reference  temperature  to  take 
the  temperature  profile  resulting  from  the  boundary  layer  velocity  distribution 
and  sample  surface  temperature  into  account. 

Mach  number  of  flow  through  the  test  sample: 


where  ps  is  the  sonic  static  pressure  if  M  =  1  or  the  downstream  pressure,  P2 
if  the  flow  is  subsonic. 

Free  stream  static  temperature: 

v  -  1  2-1 

Tg  =  Tx(l  +  ("  2~~) M  ) 

Free  stream  velocity: 

Ue  =  [5  x  104  YCV(T1  -  TE)]1/2 


Reference  temperature: 

Tre£  =  (Te  *  W/2  *  4'4  x  10 


where  T  is  the  test  sample  surface  temperature. 
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Viscosity: 

U  =  7.0  x  10’7  Tref1*5  (Tref  +  198) 


Density: 


pref  ’  PS  (BPS  *  RTref) 


-1 


Turbulent  flat  plate  heat  flux  at  the  location  of  the  in-wall  thermocouple: 


q  =  0.0S2(prefue)°-V-2  v  (Tj  -  TSAMP) 


The  sample  surface  temperature  is  determined  from  the  one-dimensional 
unsteady  state  heat  conduction  equation, 


3T 
.3 1 


•? 

3"T 


a- 


3  2 


x 


with  the  surface  boundary  condition. 


Q  = 


for  X  =  0,  where  Q  is  the  heat  flux  to  the  sample  surface.  This  equation  does 
not  consider  the  effects  of  cylindrical  geometry. 


A  finite  difference  technique  using  a  geometrical  node  grid  spacing 
was  incorporated  into  the  Calspan  code  to  solve  the  unsteady  heat  conduction 
equation.  The  general  finite  difference  relationship  is  given  by 


AT. 

i 

At 


2ct 

(F+DAx2.^ 


(T 


i-1 


T.  +  —  T  ) 
F  i  F  i+lj 


where  a  is  the  thermal  diffusivity 
T  is  the  temperature  rise 
t  is  the  time  interval 
Axj__]_  is  the  thickness  of  the  i-lst  grid 

F  is  the  geometrical  multiplier  with  the  thickness  of  the  ith 
grid  being  F  times  that  of  the  i-lst  grid. 


The  exposed  surface  boundary  condition  is  satisfied  by  first  determining  a 
fictitious  temperature  in  free  space, 
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2q 


T 

o 


Ax, 


+  T„ 


This  temperature  is  then  used  to  establish  the  surface  temperature  rise  by 


5 

At 


-A-  (T  -  2T.  +  T_) 
Ax,  2  v  o  1  2 


Barrel  Flow 


The  flow  through  the  barrel  is  expressed  as  an  input  of  mass  and 
energy  to  the  volume  between  the  test  sample  and  the  projectile.  As  mass 
and  energy  are  accumulated,  this  is  expressed  in  terms  of  pressure  and 
temperature,  which  in  turn  provides  the  accelerating  force  for  the  projectile. 

In  this  calculation,  the  quasi-steady  assumption  of  the  previous  calculation 
is  relaxed  by  allowing  pressure  acting  on  the  base  of  the  projectile  to  be 
modified  according  to  the  Mach  number  of  the  flow  at  the  base  of  the  projectile. 
In  this  way,  the  unsteady  expansion  effects  of  the  flow  through  the  barrel  is 
taken  into  account.  The  equation  of  state  and  the  basic  energy  and  mass 
conservation  equations  are  the  same  as  for  the  test  gas'  in  the  driven  tube. 

These  equations  are  used  to  define  "2"  conditions  in  the  barrel. 


Projectile  motion  is  calculated  by  an  approximate  technique  that 
involves  determination  of  the  flow  Mach  number  at  the  projectile  base. 

Specific  enthalpy: 


H 


2 


E 


2  Y 


M 


O 

C? 


where  Mo-  is  the  mass  of  gas  contained  in  the  barrel. 

c> 

Static  enthalpy  at  the  projectile  base: 

V  2 

H?  =  2.5  x  10  H7  - 

“oo  “  ~ 


where  is  the  projectile  velocity 
Mach  number: 

M  =  V  { (y  -  1)  H2  }"1/2 

P  00 

The  assumption  is  made  that  the  pressure,  p9,  is  not  the  total  pressure,  but 
is  a  static  pressure  at  the  barrel  origin. 
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It  is  further  assumed  that  the  difference  in  pressure  between  this  location 
and  the  projectile  base  is  equal  to  the  difference  between  the  pressure  and 
the  total  pressure.  This  is  expressed  as  a  Mach  number  function: 


Y-l  9  -Y/Y-1 

(1+2  M  ) 


This  allows  the  static  pressure  at  the  projectile  base  to  be  expressed  in 
terms  of  the  barrel  origin  pressure,  p^. 

From  this  the  projectile  acceleration: 


a  =  32.2  (p_  -  p  )/W 

pr  ^2ot  p 

where  p^  is  the  projectile  resistance  to  motion,  and  IV  is  the  projectile 

weight,  with  p  =  p  +  p  .  v  ” 

3  rr  rr  *r  X 

o  x  pr 

Projectile  velocity  change: 


AV  =  a  At 


pr  pr 

Projectile  displacement: 
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V 

pr 


At  +  a 


pr 


At" 

2 
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LIST  OF  CONSTANTS  AND  VARIABLES  USED  IN  THIS  SHOCK  TUBE  GUN  PROGRAM 
DRIVER  GAS  PARAMETERS 

PCH  INITIAL  DRIVER  GAS  PRESSURE  -  MPA 

TO  INITIAL  DRIVER  GAS. TEMPERATURE  -  K 

VOLO  INITIAL  DRIVER  GAS  VOLUME  -  M**3 

RDVR  DRIVER  GAS  CONSTANT  -  J/KG-K 

CVDVR  DRIVER  GAS  SPECIFIC  HEAT  -  J/KG-K 

PISTON  PARAMETERS 

AT  PISTON  FACE  AREA  -  M**2 

KP1  PISTON  DRAG  RATE  CONSTANT  -  N/M 

VELP  INITIAL  PISTON  VELOCITY  -  M/SEC 

WSHELL  PISTON  MASS  -  KG 

SHELRV  PISTON  MASS  DURING  RETURN  -  KG 

XPMAX  MAXIMUM  PISTON  DISPLACEMENT  -  M 

TEST  GAS  PARAMETERS 

P 20  INITIAL  TEST  GAS  PRESSURE  -  MPA 

T20  INITIAL  TEST  GAS  TEMPERATURE  -  K 

R2  TEST  GAS  CONSTANT  -  J/KG-K 

CV2  STANDARD  TEST  GAS  SPECIFIC  HEAT  -  J/KG-K 

CVT2  COMPUTED  TEST  GAS  SPECIFIC  HEAT  -  J/KG-K 

CVEXP  EXPONENT  FOR  COMPUTING  SPECIFIC  HEAT  OF  TESTGAS  AT  T1 

ALFA2  VAN  DER  WAALS  CONSTANT  -  MPA- M* *3 /( KG -MOL )**2 

BETA2  VAN  DER  WAALS  CONSTANT  -  M**3/KG-MOL 

GAM2  TEST  GAS  RATIO  OF  SPECIFIC  HEATS 

TEST  VOLUME  PARAMETERS 

A2  TEST  VOLUME  AREA  -  M**2 

CD2  TEST  VOLUME  EXIT  FLOW  COEFFICIENT 

DOR  1 2  TEST  SAMPLE  FLOW  CHANNEL  DIAMETER  -  M 

VOL  1 F  TEST  GAS  COLLECTION  CHAMBER  VOLUME  -  M**3 

VOL 20  INITIAL  TEST  GAS  VOLUME  -  M**3 

VOL2F  FINAL  TEST  GAS  VOLUME  -  M**3 

PROJECTILE  PARAMETERS 
BARL  BARREL  LENGTH  -  M 

BORED  BORE  DIAMETER  -  M 

PSTART  PRESSURE  ABOVE  WHICH  PROJECTILE  IS  ALLOWED  TO  MOVE  -  MPA 

RESO  Y-INTERCEPT  OF  PROJECTILE  RESISTANCE  FUNCTION  -  N 

RESS  SLOPE  OF  PROJECTILE  RESISTANCE  FUNCTION  -  N/M 

RESC  STEADY-STATE  PROJECTILE  RESISTANCE  -  N 

VPROJ  INITIAL  PROJECTILE  VELOCITY  -  M/SEC 

XPROJ  INITIAL  PROJECTILE  DISPLACEMENT  -  M 

WPROJ  PROJECTILE  MASS  -  KG 

TIME  PARAMETERS 

DELT  TIME  INCREMENT  -  SEC 

PRCI  DELT  MULTIPLIER  FOR  PRINT  INTERVAL 

TFO  LIMITING  TIME  -  SEC 
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100  FORMAT  (8F10.5) 

101  FORMAT  (ID 

102  FORMAT  ( F5 . 2 , 5A3 , 6E 1 0 . 3 ) 

103  FORMAT  <  12H1  INPUT  DATA , T48 TEST  N0.',T57,I1> 

104  FORMAT  ( ////T6 TEST  GAS  MIXTURE  BY  MOLE.  FRACTION '/ ) 

105  FORMAT  ( /T6 , F6 . 3 , T1 7 , 5A3  ) 

110  FORMAT  <///T6, ' INITIAL  CONDITIONS  FOR  TEST') 

120  FORMAT  <///T7,'PCH  = ' ,T1 6 , F 1 2 . 5 , T29 , ' PSF * , T48 , ' P20  = 1 ,T57 , F 1 2 . 5 , 
+T70 , 'PSF' , T90 , 'AT  =',T99,F12-5,T112, 'FT**2'//T7, 'TO  -',T16,F12.5, 
+T29 , ' R ' ,T48 , 'T20  =' , T57 , F 1 2 . 5 , T70 , 'R’,T90,’KP1  =',T99,F12.5,T112, 
+'LBF/FT'//T7, 'VOLO  = ' , T1 6 , F 1 2 . 5 , T29 , ' FT**3 ' , T48 , ' R2  = * , T57 , F 12 . 5 , 
+T70 , ' FT-LBF/LBM-R ' ,T90, 'VELP  »',T99,F12.5,T112, ' FT/SEC ' //T7, 

+ ' RDVR  =' , T16 , F 1 2 . 5 , T29 , 'FT-LBF/LBM-R' ,T48, ' CV2  =’ , T57 , F 1 2 . 5 , T70 , 

+  ' BTU/L  BM-R ' , T90 , ' WSHELL  =' , T99 , F 1 2 . 5 , T 1 1 2 , 'LBM'//T7, 'CVDVR  =' ,T16, 
+F12.5.T29, 'BTU/L BM-R' ,T48, ' CVT2  = ' , T57 , F 1 2 . 5 , T70 , 'BTU/L BM-R' ,T90, 

•  + 1  SHE  L  RV  = ’ ,T99 ,F 12 . 5 ,T1 12 , 1 LBM' //T48 , ' CVEXP  = '  , T57 , F 1 2 . 5 , T90 , 

+ ' XPMAX  =',T99,F12.5,T112,'FT'//T7, ' BARL  = ' , T1 6 , F 1 2 . 5 , T29 , ' FT ' , T48 , 
+ ' ALFA2  =' , T57 , F 1 2 . 5 , T70 , 'PSF-FT**3/( LB -MOL )**2 ' //T7 ,' BORED  =',T16, 
+F12.5.T29, ' FT ' , T48 , ' BETA2  =' , T57 , F 1 2 . 5 , T70 , ' FT**3/LB-MOL ' ,T90, 

+'A2  =' ,T99;F12.5,T112, 'FT**2'//T7, 'P START  =‘,T16,F12.5,T29,'PSF', 
+T48 , ' GAM2  = ' , T57 , F 1 2 . 5 , T90 , ' CD2  = ' , T99 , F 1 2 . 5//T7 , ' RE  SO  =',T16, 

+  F12.5.T29',  '  L  BF  1  ,T90,  'DOR  12  ■  '  ,T99,F  12 . 5  ,T1 12  ,  'FT'//T7,  'RESS  ='  , 

+  T 1 6 , F 1 2 . 5 , T29 , 1 L  BF /FT ' ,T4  8,'DELT  =',T57,F12.5,T70,'SEC',T90, 

+  1  VOL  1 F  ='  , T99 , F 12 . 5, T1 1 2 ,  1 FT**3'//T7,  'RE  SC  ■= '  ,  T1 6 ,  F  1 2 . 5 ,  T29  ,  '  LBF  '  , 
+  T48 , 'PRC  I  =' , T57 , F 1 2 . 5 , T90 , ‘VOL  20  ='  , T99 , F 1 2 . 5 , T 1 1 2 , ’FT**3'//T7, 

+  ' WP ROO  =',T16,F12.5,T29,'LBM',T48,'TFO  “' , T57 , F 1 2 . 5 , T70 , 'SEC'  ,T90, 
+ ' VOL2F  =',T99,F12.5,T112,'FT**3'/////> 

130  FORMAT  (//T7,'PCH  = ' , T1 6 , F 1 2 . 5 , T29 , ' MPA ' , T48 , ' P20  = ' ,T57 , F 12 . 5 , 
+T70 , 'MPA' ,T90, 'AT  = ’ , T99 , F 1 2 . 5 , T 1 1 2 , 'M**3'//T7, 'TO  = ' ,T1 6 , F 12 • 5 , 
+T29 , ' K 1 ,T48 , 'T20  =',T57,F12.5,T70,'K',T90,'KP1  =',T99,F12.5,T112, 
+'N/M'//T7, 'VOLO  =',Ti6,F12.5,T29,'M**3',T48,'R2  =',T57,F12.5,T70, 

+ ' J/KG-K ' ,T90 , ' VELP  = ' , T99 , F 1 2 . 5 , T1 1 2 , ' M/SEC ' //T7 , ' RDVR  =',T16, 
+F12.5.T29, ’ J/KG-K' ,T48, ' CV2  =' , T57 , F 1 2 . 5 , T70 , ' J/KG-K ', T90 , 

+ 1 WSHELL  =' ,T99,F12.5,T112, ' KG ' //T7 , 'CVDVR  = ' ,T16 , F 12 . 5 ,T29 , 

+ 'J/KG-K' ,T48, ' CVT2  =', T57, F12. 5, T70, 'J/KG-K', T90,'SHELRV  =',T99, 
+F12.5.T112, ‘ KG ' //T48 , 'CVEXP  = ' , T57 , F 1 2 . 5 , T90 , 'XPMAX  = 1 ,T99 , F 12 . 5 , 
+T1 12, ' M ' //T7, 'BARL  = ' , T1 6 , F 1 2 . 5 , T29 , ’ M ' , T48 , ' AL FA2  =',T57,F12.5, 

+  T70 , ' MPA-M**3/( KG-MOL  )**2  '  /  /T7 ,  ' BORED  =' ,T16,F12.5,T29, ' M '  ,T48, 

+ ' BETA2  =' ,T57,F12.5,T70, ' M**3/KG-MOL ' ,T90,'A2  ■ ' ,T99 , F 1 2 .5 ,T1 12 , 
+'M**2'//T7, 'P START  =' , T1 6 , F 1 2 . 5 , T29 , 'MPA' ,T48, ' GAM2  =' ,T57,F12.5, 
+T90 , ' CD2  =' .T99.F12.5//T7, 'RESO  = ' , T1 6 , F 1 2 . 5 , T29 , ' N ' , T90 , 

+  ' DOR  1 2  ='  , T99 , F 1 2 . 5 , T 1 1 2 ,  'M'//T7,  'RESS  =I,T16,F12.5,T29,‘N/M',T48, 
+  ' DELT  =  ' ,T57,F12.5,T70, 'SEC' ,T90,' VOL  IF  =' ,T99,F12.5,T112, ' M**3 1 // 
+T7 , ' RE  SC  ='  ,T16,F12.5,T29, ' N '  ,T48,  'PRC  I  = '  , T57  ,  F 1 2 . 5 , T90 , 

+ ' VOL20  =' , T99 , F 1 2 . 5 , T 1 1 2 , 'M**3'//T7, ' WPROJ  =' ,T16,F12.5,T29, 'KG' , 
+T48 , ' TFO  =' ,T57,F12.5,T70, 'SEC' ,T90, 'VOL2F  =' , T99 , F 1 2 . 5 , T 1 1 2 , 

+ ' M**3 '/////> 

187  FORMAT  (18H1  TABULATED  OUTPUT) 

188  FORMAT ( ///T6, 'TIME' ,T20, 'PI ' ,T34, ' T 1 ' ,T48, 'Ml ’ ,T62,  VOL  1 ' ,T76,  VEL 
+  P '  , T90 ,  ' VPROJ 1 , T104 , ‘WORK '  ,T118,'QFLUX'/T6,'PCH'  ,T20,  ' P  2 '  ,T34,  ' T2 ’ 
+  ,T48 , ' M2 '  , T62 , ' VOL  2 ' ,T76, ' XP '  ,T90, ‘ XPROJ '  ,T104, ‘ UE ' ,T118, 'QrOT'/ 
+T4 , ' TSAMP ( 1 > ' , T 1 8 , ' TSAMP ( 2 ) ' ,T32 , 'TSAMP( 3 ) ' ,T46 , 'TSAMP( 4 ) ' , 

+  T60 , ' TSAMP ( 5  )  '  ,T74,  'TSAMP(6)'  ,T88,  'TSAMP(7)‘ ,T102,  'TSAMP<8)' , 

+  T116,  ' TSAMP  <  9 ) ' // ) 

190  F ORMAT (9E14.6/9E14.6/9E14.6/ > 
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READ  <5,101  )  NUMBER 
DO  99  NM»1, NUMBER 

DIMENSION  RG( 9  > , CVG<  9 ) ,  CVTG( 9  > , CVEXPG(  9 ) ,ALFAG( 9  > , BETAG<  9  )  ,  XG( 9  ) 
DIMENSION  G1(9),G2(9),G3(9),G4(9),G5(9) 

DIMENSION  TNEW< 40 > ,TSAMP<  40  > , DELX<  40 ) , DXNEW( 40  > , DTEMP ( 40 > 

DIMENSION  TSAMPM( 40  ) 

REAL*4  KP 1 ,M1 ,M2 ,M20 
REALM  KP 1M , Ml M , M2M , M20 

DATA  VOL  0/3 1 . 5/ , CVDVR/Q . 1 76/ , RDVR/55 . 2/ 

DATA  AT /0. 3068/, WS HELL/150. 0/,SHELRV/150.0/,XP MAX/81 .0/ 

DATA  A2/0.3068/ ,CD2/0 ,75/,D0R12/0. 5/,VOLlF/0 .082/ ,VOL2F/0 .00223/ 
DATA  BARL/15.0/, BORED/O. 0984/, RESO/7200./,RESS/4320./,RESC/7200./ 
DATA  WPROO/O.25/ 

00  =  1 
KK=  1 
NN  =  0 

QFLUX=0 . 0 
UE  =  0 . 0 
R2  =  0 . 0 
CV2=0 . 0 
CVT2=0 . 0 
CVE  XP  =  0 . 0 
AL  FA2  =  0 . 0 
BETA2=0 . 0 
GAM2= 1 .0 
P  20= 1 4 . 7 
T20=530 . 0 
XP ROO  =  0 . 0 
VEL  P  =  0 . 0 
VP  ROO  =  0 . 0 
VOL 20  =  0 . 0 
DOR  1 2  =  0.5 

READ  (5,100)  PCH, TO, KP 1 , P START , DE LT , P RC I ,TFO, FACTOR 
READ  (5,100)  SIUIN,SIUOUT 
IF  (SIUIN.EQ.O.  )  GO  TO  1 
PCH  =  PCH  *  145.04 

1  READ  (5,101)  NG 
WRITE  (6,103)  NM 
WRITE  (6,104) 

DO  3  I D= 1 , NG 

READ  (5,102)  XG( ID  )  ,G1( ID ) ,G2( ID ) ,G3( ID  > ,G4( ID  >,G5(  TD ) ,RG( ID ) ,CVG( 
+ I D ) , C VTG ( ID ) ,CVEXPG( ID ) ,ALFAG( ID), BETAG ( ID ) 

IF  (SIUIN.EQ.O.  )  GO  TO  2 
RG( ID  )  =  RG( ID  )  *  0.186 
CVG(ID)  «  CVG(ID)  *  0.000239 
CVTG(ID)  =  CVTG(ID)  *  0.000239 
ALFAG(ID)  =  ALFAG(ID)  *  153746.04 
BETAG (ID)  =  BETAG(ID)  *  16.214285 

2  WRITE  (6,105)  XG< ID),G1(ID),G2( ID),G3(ID),G4( ID),G5( ID) 

R2  =  R2  +  XG ( ID  ) *RG ( ID  ) 

CV2  =  CV2  +  XG( ID  )*CVG( ID ) 

GAM 2  =  GAM2  +  XG ( I D >*< RG( I D >/( 778 . WCVG( I D > ) > 

ALFA2  =  ALFA2  +  XG( ID  )*ALFAG( ID ) 

BETA2  =  BETA2  +  XG( I D  )*BETAG( ID ) 

3  CONTINUE 

ABORE  =  0.785*  BORED**2 

CV=CVDVR 

DELTO=DELT 
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DOR12=DOR12/12.0 
A12  =  0.7854*0 OR  12**2 
I PRC= I F I X<  PRC  I  ) 

I P  R I NT  =  0 
PCH=PCH* 144. 

PCHP  =  PCH/ 1 44 . 

P  20  =  P  20* 144.0 
PARTIM  =  0. 

PSTART=PSTART*144. 

QTOT=0 . 0 

RGAS=RDVR 

CP  =  RGAS/778 . +CV 

GAM=CP/CV 

GM  =  PCH/(RGAS*TO)*VOLO 
E  =  GM* ( CP -RGAS/778  * )*T0 
TIME  =  0. 

VMAX=XPMAX*AT 
5<P  =  XPMAX 

C  INITIAL  CONDITIONS  -  UPSTREAM  -  FROM  PISTON  FACE  TO  TEST  SAMPLE 
P1-P20 

V0L1=V0L1F+XP*A2 
T1 «T20 

Ml «P 1 *VOL 1 /R2/T1 
CV I  2“0 . 0 
DO  4,J«lfWG 

CVT2«CVT2+XG<  0  )  ,VCVTG<  J  )*(  T1  -460 . 0  >**CVEXPG<  0  ) 

4  CONTINUE 

E 1 UM1 *( CV2+CVT2 )*T1 
UNO  1 -M 1 /VOL  1 
VV1 « 1 . 0/RH01 

C  INITIAL  CONDITIONS  -  DOWNSTREAM  -  FROM  TEST  SAMPLE  TO  PROJECTILE  BASE 

P  2 53  P  2  0 

V  0  L  2  «  V  0  L  2  K  +  X  P  R  0  J  *  A  B  0  R  E 
T2--T20 

M2-P  2*  VOL  2  /R2/T2 
CVT2U0 . 0 
DU  D 

»JVT2wCVT2vXG€  J  )»CVTG(  J  )*( T2-4G0 . 0 >**CVEXPG<  J  ) 

5  CONTINUE 

E2«M2vV(  CV2+CVT2  )*T2 
10102 53  M2 /VOL  2 

V  V  2  ~  1  .  /  Ill  102 
\'i\\  i  T  t;  (6,110) 

IF  (SIUUUT.KQ.O. )  GO  TO  6 
PC MM  *  PWI  *  0.00004 7 SB 
P20M  «  P20  *  0.00004780 
ATM  ~  AT  *  0.093 
TOM  “  TO  *  b . / 9 . 

T20I1  =  Y20  *  5  „  /  9  . 

KPlI'l  «  KP I  *  14. 690 
VOLOM  ^  VOLO  *  0.020 
ROM  “  R2  ft  b.376 
VLi L P 1-1  “  VELP  ft  0.206 
RDVRI'l  “  RDVU  *  5.376 
CV2M  =  CV2  *  4134.1 
WoDlILM  «  V/3HE LL  *  0.454 
CVDVRM  «  CVDVK  *  4184.1 
CVT2U  «  CVT2  *  4104. 1 
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SHELRM  =  SHELRV  *  0.454 
XPMAXM  =  XPMAX  *  0.305 
BARLM  =  BARL  *  0.305 
ALFA2M  =  ALFA2  *  0.0000065 
BOREDM  =  BORED  *  0.305 
BETA2M  =  BETA2  *  0.061674 
A2M  =  A2  *  0.093 
PSTARM  =  PSTART  *  0.00004788 
RESOM  =  RESO  *  4.45 
DOR12M  =  DOR12  *  0.305 
RESSM  =  RESS  *  14.596 
VOL  1 FM  =  VOL  1 F  *  0.028 
RESCM  =  RE  SC  *  4.45 
VOL20M  =  VOL20  *  0.028 
WP ROOM  =  WPROO  *  0.454 
VOL2FM  =  VOL2F  *  0.028 

WRITE  (6,130)  PCHM.P20M ,  ATM ,T0M,T20M,KP1M, VOLOM ,R2M,VELPM,RDVRM, 
+CV2M,WSHELM,CVDVRM,CVT2M, SHELRM, CVEXP , XPMAXM , BARLM ,ALFA2M , BOREDM , 
+  BETA2M,A2M,PSTARM,GAM2,CD2.,  RE  SOM,  DOR  12M,  RESSM,  DELT,  VOL  1FM,  RESCM, 
+PRCI , VOL20M,WPROJM,TFO, VOL2FM 
GO  TO  7 

6  WRITE  (6,120)  PCH , P20 , AT , TO , T20 , KP 1 , VOLO , R2 , VE L P , RDVR , CV2 , WSHELL , 
+CVDVR , CVT2 , SHELRV, CVEXP , XPMAX , BARL , ALFA2 , BORED , BETA2 , A2 , PSTART , 
+GAM2, CD2, RESO, D0R12, RESS, DELT, VOL  IF , RE  SC, PRC  I , VOL  20, WPROO ,TFO, 
+V0L2F 

7  DEPTH  =  0.5/12.0 
XK  =  19.3/3600. 

ALPHA  =  XK/57.6 

F  =  1 . 3 

DELXO=SQRT< ALPHA*DELT/0 .25 ) 

SUMX=DELXO 
D  E  L  X ( 1 )  =  DE  L  XO 
DO  8  1=1,40 
TSAMP ( I )=T0 
TNEW( I )=T0 

8  CONTINUE 

DO  9  1=1,39 
D E L X ( 1  +  1 >  =  D E L X ( I  )*F 
IF( SUMX.GE .DEPTH  )  GO  TO  9 
S  U  M  X  =  S  U  M  X  +  D  E  L  X  (  I  ) 

KTEMP  =  I 

9  CONTINUE 

CONST1 =ALPHA"2 .0/(1. 0  +  F  ) 

CONST2  =  ( 1 .0  +  F  >/F 
CONST3= 1 . 0/F 

TZIP=QFLUX*2 . 0*DELX(  1  )/XI<+TSAMP(  2  ) 

10  CONTINUE 
T=T0 

RGAS=57 . 4+ . 333E -5*PCH 
CV=0. 1 75+ . 1 83E-4*T 
CP=RGAS/778 . +CV 
GAM=CP/CV 

I  F ( VELP.LT.-O. 1  )  WSHELL  =  SHELRV 

C  COMPUTE  THE  DYNAMICS  AND  THERMODYNAMICS  OF  PISTON  MOTION.  . 

24  X D  =  X P M A X - X P 

D  RAG  =  KP 1 *A  B  S ( V  E  L  P  ) 

FOVM  =  (  (  PCH ••AT-P  1  ftA2  > -DRAG  ) /WSHELL *32 . 17 
26  DVELP  =  FOVM  *  DELT 
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28  VELP  =VE L P  +DVELP 

I F ( XP . LE . 0 . 0000 1 . AND. VELP .GT. 0.0  )  VELP  =  0.0 
DXP=DELT*(VELP+F0VM*DELT/2.) 

XP  =XP -DXP 

I F ( XP  .  L E  . .00001  )  GO  TO  99 
WORK  =  PCH  *  AT  *  DXP/778. 

E  =  E  -  WORK 

30  VOL  =  VOL 0+AT*XD 
RHOCH  =  GM/VOL 

T  =  E/CCCP-RGA S/778.  )*GM> 

PCH  =  RHOCH*RGAS*T 
C  COMPUTE  TEST  GAS  PROPERTIES. 

CVTX=0 . 0 
CVTY-0 . 0 
DO  31  0=1, NG 

CVTX=CVTX+XG< J  )*CVTG( J  >*<T2-460.0  >**CVEXPG< J ) 
CVTY=CVTY  +  XG( J  >*CVTG< J  )*(Tl-460 .0  )**CVEXPG( J ) 

31  CONTINUE 

CVY= ( (CV2+CVTY)*T1-ALFA2*RH01/778.0)/T1 
CVX  =  < ( CV2  +  CVTX  )*T2-ALFA2*RH02/778 .0  >/T2 
RESIST  =  RESO  +  RESS  *  XPROJ 

32  IF( RESIST.LT. RE SC)  RESIST=RESC 
H2«E2*GAM2/M2 

H  2  I N  F  =  2  5  0  0  0 . *  H  2 - V  P  RO J  *  *  2 / 2  . 

XM2INFBVPROJ/SQRT<  (GAM2-1  .0>*H2INF  ) 

F  M  =  (  1 . 0  +  ( GAM 2- 1 . 0  )/2 . 0*XM2 I NF**2  )**( -GAM2/<  GAM2- 1 .0 ) ) 
DELP«P2*<  1 . 0-FM  >/<  1  ♦  0  +  F M  ) 

33  P20=:P2  +  DELP 
P2INF=P2-DELP 

PFORCE  =  ABORE  *  ( P2INF-RESIST) 

IFCPFORCE.LT. 0.0)  P FORCE  =  0 . 0 

IFCP2.LT. P START  )  PFORCE=0.0 

APROO  «  PKORCE/WPROJ*32. 17 

DXPROJ  =  VPROO  *  DELT  +  APROO  *  DELT*  *2/2 . 

35  XPROO  «  XPROJ  +  DXPROJ 

IFCXPROO.GE.BARL)  GO  TO  99 
V  PRO  Cl  «  VP  110  J  +  APROO  *  DELT 
VOL  1« VOL  IF *XP *5fA2 
VOL  2“ VOL  2F  +XPROO *AUORE 
3G  P1AVE-P1 
X  Ml « 1-11 
XM2«M2 
I  P  A  S  S  =  0 
P 1 SAV-P 1 
P  2SAV-P2 
37  CONTINUE 

I  PASS- 1  PASS* 1 

P  GTAU-P  1 V/  (  1  .0+CGAM2-1  .0  )/2.0  )**< -GAM2/C GAM2-1 .0)  ) 
TTOT-T1 
P  S“P  STAR 
P  TOT = P 1 

1 1:  (  P 2 0  .  G T  .  P  1  )  GO  TO  47 
IFC  PSTAR.GE . P20 )  GO  TO  38 
PS“P20 

D  E L  TM 1 “ 0 . 0  2  * A 1 2  * SORT  C C P1-P20 >*P 1/R2/T1 )*DELT 
GO  TO  49 

3G  TSTARbT1 *2 • 0/< GAM2  + 1.0) 

DELTM1UPSTAR*A12V,<JQRT< GAM2tt32 . 2/R2/TSTAR  >*DELT 
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GO  TO  49 

47  CONTINUE 
PSTAR=PSTAR*P20/P1 
TTOT =T2 
PS=PSTAR 
PTOT=P20 

I F ( P STAR . GE . P I )  GO  TO  48 
PS  =  P1 

DELTM1=-8.02*A12*SQRT( (P20-P1 )*P20/R2/T2 )*DELT 
GO  TO  49 

48  TSTAR=T2*2.0/<GAM2+1 .0) 

DELTM1=-PSTAR*A12*SQRT( GAM2*32 .2/R2/TSTAR)*DELT 

49  CONTINUE 
XM1=M1-DELTM1 

■  XM2=M2+DELTM1 
RH01=XM1/V0L 1 
RH02  =  XM2/V0L  2 
VV1=1 .0/RH01 
VV2=1./RH02 

WORK  =  -ABORE"DXP  ROO*P  2  INF/778. 0-P2/288. *XPROO  *3 . 1416*BORED*DELT 

IFIVPROO.LT.IOO. )  GO  TO  58 

DELH=CVY*DELTM1*T1*GAM2 

50  I F ( DELTM1 . LT. 0 . 0 )  DELH=  CVX*DELTM1*T2*GAM2 
E  X  =  E  2+V/ORK+DE  L  H 

T2=EX/CVX/XM2 

PX  =  <  R2*T2  )/<  VV2-BETA2  )+ALFA2/VV2J,I*,2 
P  2 “ < PX+P2SAV )/2 . 0 
E Y=E 1 -DELH+P 1AVE/778 . 0*A2*DXP 
T1=EV/CVY/XM1 

PY  =  ( R2*T1  >/( VV1-BETA2 )+ALFA2/VVl**2 
P  1  -  (  P  V + P 1 S AV  > / 2 . 0 

54  DELP1-PV-P1 
P1AVE=P1+DELP 1/2.0 
P20-P2  +  DEL.  P 

55  I  F  < IPASS.LT. 3)  CO  TO  37 
E  1  =  E  V 

E2  =  EX 
P  1  =  P  V 
P  2  =  P  X 
M1=XM1 


M2=XM2 

XMOR=SQRT( ( ( PTOT/PS )**( (GAM2-1 .0  >/GAM2  >-l . 0 )*2 . 0/< GAM2-1 . 0  >  > 
TE=TTOT /( 1 . 0+( GAM2-1 . 0 )/2 . 0*XMOR**2  ) 

UE  =  SQRT (  5 . 0E4*CVY*GAM2"( TTOT-TE  )  ) 

TRE F  =  (  TE+TSAMP <  1  )  )/2 . 0+4 . 4E-6*UE**2/CVY/GAM2 
VIS=7.0E-7*TREF**1 . 5/< TRE F  +  1 98 . 0  > 

RFIOREF  =  PS/(BETA2*PS  +  R2*TREF  ) 

QF L UX  =  0 . 052*(  RHOREF*UE  )<,*0. 8*VIS**0 . 2*CVY*GAM2*(  TTOT-TSAMP  (  1  )  > 
I F ( TSAMP ( 1 >.GT. 3100.0)  GFLUX=QFLUX* < TTOT-3100 .0 )/( TTOT-TSAMP ( 1 
QFLUX  =  QFLUX"F ACTOR 


TZIP=QFLUX*2 . 0*DELX( 1 > /XK+TSAMP ( 2  ) 

DTEMP ( 1  )=ALPI!A/( D E L X ( 1  )**2  ) * <  1 2 1 P - 2 . 0*TSAMP( 1 )+TSAMP<  2  )  > 

DO  56  K-2 , KTEMP 

DTEMP ( K I^CONSTl / ( DELX { K- 1 >**2 )*(TSAMP ( K-l )-CONST2*TSAMP(-K )  + 
*  CON ST3* TSAMP ( K+ 1 ) ) 


)  > 


56  CONTINUE 

DO  57  I<=1,  KTEMP 

TSAMP ( K)=TSAMP<  IO+DTEMP< K)*DELT 
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57  CONTINUE 

TSAMP ( KTEMP  +  1 )=TSAMP( KTEMP-1 ) 

QT0T=3QT0T  +  QFLUX*DELT 
IF( P20.GT.P1  )  UE  =  -UE . 

GO  TO  60 

58  CONTINUE 

VOL  =VOL 1 +VOL  2 
XMTOT  =  M 1 +  M2 

ESUM=El+E2+Pl*A2*DXP/778.0  +WORK 

Ml aXMTOT* VOL  1 /VOL 

M2=XMTOT*VOL2/VOL 

El=ESUM*VOLl/VOL 

E  2  =  E  SUM* VOL  2/ VOL 

Tl=ESUM/CVX/XMTOT 

T2  =  T  1 

VVT= VOL / XMTOT 

Pl  =  { R2*T1  )/( VVT-BETA2 ) +AL FA2/VVT**2 
P  2  =  P  1 

RH01=M1/V0L1 
RH02  =  M2/ VOL  2 
VV1“1 . 0/RH01 
VV2= 1 . 0/RH02 
60  CONTINUE 

IF  ( VOL  1 ,  GT . 0 . 2*VMAX  )  GO  TO  80 

I F ( VOL  1 • LE • 0 • 2*VMAX ) DELT®0 . 1 *DELTO 

I K  <  VOL  1 .  LE  .  0 . 2*VMAX  .  AND  .  1<I<  .  EQ . 1 )  GO  TO  62 

I F ( VOL  1 .  L  E  .  0 . 0  2  * VMAX )DELT=0.01*DELT0 

I F ( VOL  1 . L  E ♦ 0 . 0  2 * VMA X . AN D . KK ♦ E Q . 2 )  GO  TO  62 

GO  TO  80 

C  REVISE  GRID  SIZE  DUE  TO  CHANGE  IN  TIME  STEP. 

62  DELXO-SQRT <ALPHA*DELT/0.25  ) 

SUMX-DE  L  XO 

DXNEW<  1  )“DELX0 

DO  64  1=2,40 

DXNEW< I }°DXNEM< 1-1 >*F 

IF<8UMX.GE. DEPTH)  GO  TO  64 

SUM XuSUMX*DXNEl/(  I  ) 

KNMAX  «  I 
64  CONTINUE 

XKOL  D'-DELX  (  1  ) 

XKNEW=0 . 0 
K»  1 

KNEW® 1 

A={ TZIP -TSAMP <  2 ) )/( 2 . 0*DELX( « ) ) 

B “ (  TSAMP  (  1  ) -TSAMP  (  2  ) -A*DE L i( (  «  ) ) /DEL X (  K  >* *2 
THEW< 1  )  -  TSAMP ( 1  ) 

66  KNEW«KNEW*1 

XKNEWbXI(NEW*DXNEW( KNEW- 1 ) 

TNEW(  KNEW) “TSAMP (  1  ) - (  A*XKNEW*B*XKNEW**2  ) 

1F< XKNEW.LT. XKOLD)CO  to  66 
6  7  K“l\+ 1 

XKOL l)*XKQLD*Dt:LX ( l< ) 

A«  (  TSAMP  (  K- 1  )  -TSAMP  ( !<+  1  )  >  /  (  (  F *  1 . 0  )  * D E  L  X  (  K-  1  >  > 
B u  (  TS AM P  <  K  >  “TSAMP  (  K+  1  > -AftD E L  X <  l<  >  >  /D E L X  <  K  > * * 2 
68  KNEW® KNEW* 1 

XKNEU-XKNEW+DXNEl/<  KNGW-l  ) 

KDEL«XKNEW-XKOLD  +  DEL>((IO 

TNEU<  KNEW  >  “TSAMP  (  K  >-< A*XDEL+B*XDEL**2  ) 
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IF( XKNEW.LT. XKOLD. AND. KNEW+1 .LE .KNMAX)  GO  TO  68 
IF( K+l .LE . KTEMP .AND. KNEW+1 .LE . KNMAX >  GO  TO  67 
DO  70  1=1, KNMAX 
DELX( I >=DXNEW( I > 

TSAMP ( I ) =TNEW( I ) 

TNEW( I )=T0 
70  CONTINUE 
KTEMP=KNMAX 
KK  =  KK+1 
80  CONTINUE 

C  PRINT  OUT  COMPUTED  RESULTS. 

IF  ( J J . EQ . 1  )  GO  TO  88 

87  NN= 1 

PARTIM=PARTIM+DELT 
I F ( TIME .GT.TFO  )  GO  TO  99 
I F  < IPRINT.LT. IPRC)  GO  TO  95 
GO  TO  89 

88  WRITE  (6,187) 

WRITE  (6,188) 

JJ=JJ+1 

89  PCHP=PCH/144. 

P1P  =  P  1/1  44.0 
P2P=P2/144 . 

T I  ME  =T I  ME  +  PART  I M 
PART  I M  =  0. 

IF  (SIUOUT.EQ.O.  )  GO  TO  92 
P1PM  =  PIP  *  0.0058946 
TIM  =  T1  *  5./9. 

M1M  =  Ml  *  0.454 
VOL 1M  =  VOL  1  *  0.028 
VELPM  =  VELP  *  0.305 
VPROJM  =  VPROJ  *  0.305 
WORKM  =  WORK  *  1.356 
QFLUXM  =  QFLUX  *  0.0113566 
PC II PM  =  PCHP  *  0.0068946 
P2P M  =  P2P  *  0.0068946 
T21-1  =  T2  b./9. 

M2M  =  M2  *  0.454 
VOL2M  =  VOL 2  *  0.028 
X P M  =  XP  *  0.305 
XPROJM  =  XPROJ  *  0.305 
UEM  =  UE  *  0.305 
QTOTM  =  QTOT  *  0.0113566 
DO  90  1=1,9 

TSAMPM(I)  =  TSAMP ( I )  *  5./9. 

90  CONTINUE 

WRITE  (6,190)  TIME ,P  1PM, TIM, M1M, VOL 1M, VELPM, VP ROJM, WORKM, QFLUXM, 

+  PCHPM , P2PM , T2M , M2M , VOL2M , XPM , XPROJM , UEM , QTOTM , ( TSAMP M( I ) , I -1 . 9  > 

GO  TO  93 

92  WRITE ( 6 , 190  >  TIME, PIP, Tl, Ml, VOL  1 , VELP , VPROJ .WORK, QFLUX, PCHP, P2P, T2 
+ ,M2,V0L2,XP ,XPROJ ,UE , QTOT , ( TSAMP ( I > , 1=1 ,9) 

93  IF(MN.EQ.O)  GO  TO  87 
IPRINT  =  0 

95  I P  R I NT= I P  R I NT+ 1 
GO  TO  10 
99  CONTINUE 
STOP 
END 
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'  '  ^  sample  printout  of  the  STG  program  follows.  Definitions  of 

"Initial  Condition"  variables  are  given  in  the  preceding  program  printout. 
However,  some  explanation  of  the  "Tabulated  Output"  variables  is  required. 


TIME 


PCH 


V0L1 


PI,  Tl,  Ml 

V0L2 


P2,  T2 ,  M2 
VELP,  XP 


Time  after  piston  release  -  seconds. 

Time  is  initially  zero  and  increases  to  a  maximum  value 
either  when  the  projectile  exits  the  barrel  or  when  it  equals 
a  limiting  value,  i.e.,  TFO. 

Driver  gas  pressure  -  MPa. 

PCH  is  a  prescribed  maximum  at  TIME  zero  and  decreases  as 
the  driver  gas  displaces  the  unlatched  piston. 

Test  gas  volume  included  from  piston  face  to  test  specimen 
inlet  -  m3. 

V0L1  is  initially  the  entire  volume  of  the  driven  tube  but 
decreases  to  the  volume  of  the  test  gas  collection  chamber 
when  the  piston  has  been  fully  displaced. 

Test  gas  pressure,  temperature  and  mass  associated  with 
V0L1 ,  measured  in  MPa,  °K  and  kg,  respectively. 

Test  gas  volume  included  from  test  specimen  inlet  to 
projectile  base  -  m3. 

V0L2  has  a  minimum  value  of  the  test  specimen  bore  volume 
at  TIME  zero  and  increases  with  projectile  displacement  to 
include  the  entire  barrel  volume. 

Test  gas  pressure,  temperature  and  mass  associated  with  V0L2, 
measured  in  MPa,  °K  and  kg,  respectively. 

Piston  velocity  and  displacement  -  m/sec,  m. 


VPROJ,  XPROJ  Projectile  velocity  and  displacement  -  m/sec,  m. 


WORK 


UE 

QFLUX 


Work  performed  by  driver  gas  on  the  driven  piston,  and  by  test 
gas  on  the  projectile.  During  piston  rebound,  the  program  also 
computed  negative  work  done  on  the  piston  by  the  test  gas  -  J. 

Test  gas  free  stream  velocity  -  m/sec. 

2 

Test  sample  surface  heat  flux  -  J/mm  -sec. 


QT0T  Total  integrated  heat  input  to  test  sample  surface  -  J/mm2. 

TSAMP  Test  Sample  surface  temperature  -  °K. 

TSAMP(l)  is  the  computed  temperature  on  the  test  sample  surface. 
TSAMP (2)  through  TSAMP (9)  are  subsurface  temperatures  computed  at 
depths  printed  in  the  nonlinear  DELX  array,  which  is  amended  when 
the  time  increment  DELT  is  changed. 
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INPUT  DATA 


TEST  NO.  1 


TEST  GAS  MIXTURE  BY  MOLE  FRACTION 


0.455 

CARBON  MONOXIDE 

0.545 

ARGON 

INITIAL 

CONDITIONS  FOR 

TEST 

PCH  - 

4.41262 

MPA 

P20  - 

0. 10135 

MPA 

AT  - 

0.02853 

M*  *3 

TO  - 

297.22217 

K 

T20  = 

294.44434 

K 

KP1  = 

87.57599 

N/M 

VOLO  = 

0  88200 

M*  *  3 

R2  « 

248.25496 

J/KG-K 

VELP  * 

0.0 

M/SEC 

RDVR  - 

296.75513 

J/KG-K 

CV2  - 

511  .26343 

J/KG-K 

WSHELL  - 

68.09999 

KG 

CVDVR 

»  736.40186 

J/KG-K 

CVT2  = 

2.29213 

J/KG-K 

SHELRV  = 

68.09999 

KG 

CVEXP'  = 

0.0 

XPMAX  » 

24.70499 

M 

BARL  = 

4.57500 

M 

AI.FA2  ■= 

0.00470 

MPA-M**3/( KG 

-MOL )**2 

BORED 

-  0.03001 

M 

BETA2  = 

0.00107 

M**3/KG-MOL 

A2  = 

0.02853 

M**2 

PSTART 

=  0.34474 

MPA 

GAM2  = 

1 . 48575 

CD2  * 

0.75000 

RESO  = 

32039.9961 

N 

DOR  1 2  * 

0.01271 

M 

RESS  - 

63054 .7187 

N/M 

DELT  « 

0.00100 

SEC 

VOL  IF  * 

0 . 00230 

M**3 

RE  SC  = 

32039.9961 

N 

PRC  I  = 

33.00000 

VOL20  « 

0.0 

M**3 

wproj 

-  0.11350 

KG 

TFO  = 

0.25000 

SEC 

VOL2F  * 

0.00006 

M**3 

TABULATED  OUTPUT 

TIME 

TSAMP ( 1 ) 

PI 

P  2 

TSArIP  (  2  > 

T1 

T2 

T  SAMP ( 3 ) 

Ml 

M2 

TSAMP ( 4  > 

VOL! 

VOL  2 

TSAMPf 5  ) 

VELP 

XP 

TSAMP ( 6 ) 

VPROJ 
XPROJ 
TSAMP <  7  > 

WORK 

UE 

TSAMP ( 8  ) 

OFLUX 

OTOT 

TSAMP <  9  > 

0.0 

0.439334E+01 
0.29  7222E+03 

0.  1  0 1 905E+00 
0. 1 0 1  905E  +  00 
0. 297222E+03 

0.294824E+03 

0.294824F+03 

0 . 297222E  +03 

0 . 979020E+00 

0 . B75736E -04 

0 . 297222E  +  03 

0.698041E+00 

0 . 624400E-04 

0 . 297222E  +  03 

0. 180703E  +  01 
0.247023E+02 
0.297222E+03 

0.0 

0.0 

0 . 297222E  +03 

0.0 

0.0 

0.297222E+03 

0.0 

0 . 0 

0.297222E+03 

0.339999E-01 
0. 41968 1E+01 
0.297222E+03 

0  .  1  OB  7  7  3E  +  00 
0.  10U773E  +  00 

0 . 29  7222E  +  03 

0 . 30 1 027E+03 

0 . 30 1 027E  +  03 

0 . 297Z22E  +  03 

0. 97900?E  +  00 

0 . 915380E-04 
0.297222E+O3 

0 . 667796E  +  00 
0.624400E-04 

0 . 297222E  +  03 

0.589426E+02 

0 . 236284E  +  02 
0.297222E+03 

0.0 

0.0 

0 . 297222E  +  03 

0.0 

0.0 

0 . 297222E  +  03 

0.0 

0 . 0 

0.297222E+03 

0 . 669999E -0 1 

0 . 3732 1 3E  +  0 1 
0.297222E+03 

0. 13I2G8E+00 
0. 1 3 1 26QE+00 

0 .297222E+03 

0.31 9595E+03 
0.31 9595E+03 

0 . 297222E  +03 

0 . 978970E  +  00 
0.  1  04009E-03 
0.297222E+03 

0.587706E+00 

0 .624400E-04 
0.297222E+03 

0.  I08QG0E+03 
0. 207849E  +  02 
0. 297222E+03 

0.0 

0.0 

0 . Z97222E+03 

0.0 

0 . 0 

0.297222E+03 

0.0 

0 . 0 

0 , 297222E+03 

0 . 999998E -0 1 

0 . 3  1  7230E  +  0 1 

0 . 297222E  +  03 

0. 1G464IE+00 
0. 1 8464 !  E  +  00 

0 . 297222E+03 

0 . 356056E+03 

0 . 356056E+03 

0 . 29  7222E  +  03 

0 . 978927E  +  00 

0  .  131208E-03 
0. 297222E  +  03 

0.465856E+00 
0.624400E-04 
0. 297222E+03 

0,  I49166E+03 

0  .  1  64586E  +02 
0.297222E  +03 

0.0 

0.0 

0 . 297222E  +  03 

0.0 

0.0 

0 .297222E+03 

0.0 

0.0 

0.297222E+03 

0 . 1  33000E+00 

0 .264352E+01 

0 . 297222E  +  03 

0.332495E+00 
0.332495E+00 
0. 297222E+03 

0  4  28  166E  +  03 
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0,33004*6+03 

0.83002  36  +  02 
0.2546616+00 
0. 3315 3G E+03 

0. 10051 1E  +  03 
0 . 284897E  +  00 
0 . 333487E  +  03 

0.  1  23510E  +  03 
0.321 7906  +00 
0 . 335R63E  +  03 

0.  I  53830E  +03 
0. 367480E+00 
0. 33868  IE  +  03 

0. 193106E  +  03 
0 . 42466  I E  +  00 
0 . 341996E+03 

0 . 240493E  +  03 
0 . 4  96280E  +  00 
0 . 345903E  +  03 

0 . 286722E  +03 
0 . 583784E  +  00 
0 .350530E  +  03 

0 .307826E+03 
0. 6830906+00 
0 . 356083E+03 

0.280861E+03 

0.781456E+00 

0.362774E+03 

0.21  9230E  +  03 
0  .  R64  1  loE+OO 
0 . 370064E  +  03 

0. 156819E  +  03 
0. 9255286  +  00 
0 . 380559E  +  03 

0.  J09783E  +  03 
0 . 96881  7E  +  00 
0.39I915FMJ3 

0.7787436+02 
0 . 999262F.  +  00 
0. 4O477/F.  +  03 

0. 566566E+02 
0.1021 1 l F  +  0 1 
0.41879 3 E+03 

0.42361 7 E+02 
0. I03723E+01 
0 . 4  3350GE  +  03 

0 . 326524E  +  02 
0. 1  04950E  +01 

0 . 44R44'H"*03 

0.246938E+02 
0.  105800E  +  0 1 
0.4631 7 6E+03 

0.  179797E  +  02 
0.  I06505E  +  01 
0. 477384E+03 

0.  1  18285E  +  02 
0.  107074E  +  01 
0 . 490029F  +  03 
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APPENDIX  II 


ISENTROPIC  EQUILIBRIUM  COMBUSTION  CODE 


Adiabatic  Compression  Program 

An  existing  equilibrium  combustion  program  has  been  modified  to 
calculate  the  properties  of  a  gas  mixture  at  various  points  on  an  isentrope. 
The  program  uses  a  list  of  thermal  properties  of  possible  constituents  of  the 
mixture.  The  program  writes  output  containing  the  composition  and  thermal 
properties  of  the  gas  mixture  at  the  input  pressures. 

Input  to  the  program  consists  of  "REACTANT"  cards  specifying  the 
composition  at  the  initial  point. on  the  isentrope  and  a  namelist,  "SINPT2," 
which  specifies  the  temperature  and  pressure  at  the  initial  point  and  a 
schedule  of  other  pressures  (up  to  26  pressures  including  the  initial  value) . 
OMIT  cards  may  be  used  if  desired  to  exclude  certain  species  from  considera¬ 
tion  as  mixture  constituents. 

The  program  operates  as  follows : 

(1)  The  initial  mixture  as  specified  by  the  reactant  cards  and  the 
mixture  ratio  specified  by  §INPT2  is  analyzed  to  determine  which 
chemical  elements  are  present.  The  number  of  gram-atoms  of  each 
element  is  calculated  as  well  as  the  mixture  molecular  weight 
and  mixture  total  enthalpy  (Subroutine  REAC) . 

(2)  The  taped  thermal  data  is  searched  (SEARCH)  and  the  names  and 
thermal  properties  of  possible  compounds  which  could  be  formed 
from  the  available  atoms  are  extracted.  This  list  is  then  com¬ 
pared  to  the  OMIT  cards  and  the  net  list  of  species  to  be 
considered  is  printed. 

(3)  Subroutine  EQLBRM  then  varies  the  composition  seeking  to  minimize 
the  Gibbs  free  energy  of  the  mixture  subject  to  the  constraint 
that  the  numbers  of  atoms  of  the  different  elements  do  not  vary. 

In  this  process,  species  with  molefractions  of  less  than  10  ^  at  any 
iteration  are  dropped  from  consideration.  Occasionally,  the  constraint 
equations 


where  a^  is  the  number  of  atoms  of  type  i  in  species  j ,  become  linearly 
dependent'.  In  this  case,  the  program  will  print  the  message  "SINGULAR  MATRIX" 
and  then  calculate  the  thermal  properties  of  the  input  mixture  at  the  specified 
pressure  and  temperature  (Subroutine  HCALC) .  Sometimes  the  mixture  free  energy 
will  have  a  very  broad  minimum  or  will  possess  two  minima 
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as  in  the  sketch.  If  this  is  the  situation,  the  program  will  print  the 
message 

"35  ITERATIONS  DID  NOT  SATISFY  CONVERGENCE  REQUIREMENTS  FOR  THE 
POINT 


and  then  calculate  the  thermal  properties  of  the  input  mixture, 

(4)  One  way  or  the  other,  the  program  will  obtain  the  entropy  at  point  #1 , 
This  entropy  will  then  be  held  constant  (SP=.TRUE-),  EQLBRM  will 
calculate  the  equilibrium  composition  at  an  estimated  temperature 

and  the  new  pressure,  calculate  the  mixture  entropy,  and  compare  it 
to  the  entropy  at  point’ #1.  If  the  two  entropies  do  not  agree,  the 
temperature  estimate  is  changed  and  a  new  equilibrium  calculated. 

(5)  If,  at  any  point,  the  equilibrium  calculation  fails,  the  control 
program  THERMP  will  call  subroutine  FR0ZEN  to  obtain  the  thermal  . 
properties  corresponding  to  the  fixed  entropy,  the  new  pressure 
and  the  composition  at  the  previous  point. 

(6)  After  15  points  have  been  calculated,  the  program  will  print  the 
results.  Then  it  will  complete  and  output  the  remainder  of  the 
pressure  schedule. 

(7)  After  completing  the  pressure  schedule,  the  program  will  look  for 
a  new  mixture  ratio  MIX. 
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REACTANTS 

C  1.0000  0 

1 

.0000  0.0 

0.0 

00 

0.0  0.1000 

M  0.0 

G 

298. 

150  F  0. 

.0 

H  2.0000 

0 

.0  0.0 

0.0 

00 

0.0  0.0500 

M  0.0 

G 

298. 

150  F  0. 

.0 

N  2.0000 

0 

.0  0.0 

0.0 

00 

0.0  0.0500 

M  0.0 

G 

298. 

150  F  0. 

.0 

AR  1.0000 

0 

.0  0.0 

0.0 

00 

0.0  0.8000 

M  0.0 

G 

298. 

150  F  0. 

.0 

OMIT 

C(S) 

NAMELISTS 

&  I NPT2 

KASE  = 

7, 

r-  300.000000 

,  .0 

,  .0 

,  .o 

•  .o 

,  .  0 

.0 

.0 

.  .0 

.0 

.0 

.0 

.0  , 

.0 

.0 

,  .0 

, 

.0  , 

.0 

.0 

.0 

.0 

.0 

.  .0 

, 

.0 

.0 

.0 

.0 

.0 

.0 

,  .0 

t 

.0 

.0 

.0 

.0 

.0 

.0 

,  .0 

.0 

.0 

.  0 

.0 

.0 

.0 

.  .0 

.0 

.0 

.0 

.0 

.  0 

.0 

.  .0 

.0  ,P 

*  14.6955996 

,  2900.80005 

,  5801.60156 

0702.39044 

11603. 1992 

.  14504.0000 

17404.8008 

20305.6016 

23206 

.3984 

26107.1992 

29008.0000 

31908.8008 

.  34809.6016 

37710.3904 

40611 .1992 

43512 

.0000 

46412.0008 

49313.6016 

52214.3984 

,  55115.1992 

, 

50016.0000 

60916.8008 

63817 

.6016 

66718.3750 

69619.1875 

72520.0000 

,  P  S I AeT ,MMHGa 

F,NSQM=F,V-  .0 

.  -0 

.0 

,  .0 

.  0 

.0 

.  .0 

.0 

.0 

.0 

.0 

.0 

.0 

,RHO-  14.695 

5996* 

,  2900.80005 

,  5801,60156 

,  8702.39B44 

,  11603.1992 

14504.0000 

17404.8008 

,  20305.6016 

23206.3984 

26107.1992 

29008 

.0000 

31908.8008 

34809.6016 

37710.3984 

.  40611.199? 

43512.0000 

46412.8008 

49313 

.6016 

52214.3984 

55115.1992 

58016.0000 

,  60916.8008 

63817.6016 

66718.3750 

69619 

.1075 

72520.0000  , ERATIO^F ,OF  « 

T,FPCT=F ,FA  =  F 

.MIX 

=  .  0 

♦  .0 

,  .0 

.0 

,  .0 

« 

.0  , 

.0 

.0 

t  .o 

, 

.0 

.0 

.0 

, 

.0 

.  0 

.0 

,TP=*T  ,  HP~F , SP 

“F  ,  TV- 

=  F.UV-F.SV  =  F,RKT«=T 

,  SHOCK  =  F , DETN3 

f! 

OTTO=F . 

CR=  .0 

,S0^  .0 

SO=*  .0 

, IONS*F , IDEBUG 

I*  0, 

TRACE 

*  .0 

SIUNIT-F ,EUNITS= 

F,TRNSPT=F,F  ROZN  =  F 

,  PUNCH  =  F , NODATA  =  F . D I F  =  F 

4EN0 


NO  INPT2  VALUE  GIVEN  FOR  OF,  EQRAT,  FA,  OR  FPCT 


SPECIES  BEING  CONSIDERED  IN  THIS  SYSTEM 

CH20 

L  5/66 

AR 

J  3/61  C 

012/67  CH 

J 1  2/72 

CH2 

0  3/61 

0  6/69 

CH3 

0  3/61  CH4 

J  6/69  CN 

0  6/66 

CNN 

012/70 

CN? 

0  9/65 

CO 

J  9/65  C02 

J  1  2/69  C2 

J  3/67 

C2H 

0  3/61 

C2M2 

J  9/65 

C2H4  ‘ 

L  5/72  C2H6 

0  3/67  C2N 

J  3/61 

C2N2 

0  9/66 

C20 

0  12/69 

C3 

0  6/68  C30? 

012/69  C4 

012/69 

C5 

J  3/77 

H 

L 1 2/69 

HCN 

J 1 2/70  HCO 

012/70  HNCO 

J  3/63 

HNO 

J  6/63 

HN02 

J  6/63 

HN03 

J  3/64  H02 

0  3/77  H2 

LI  1/65 

H20CS) 

LI  1/65 

H?0<  L ) 

J  3/61 

H20 

L  2/69  H202 

0  3/77  N 

012/70 

NCO 

012/71 

NH 

012/65 

NH2 

0  9/65  NH3 

J  6/63  NO 

0  9/64 

N02 

J  l  2/6  4 

N03 

0  3/77 

N2 

012/65  N2H4  . 

J 1 2/6  4  N20 

0  9/64 

N204 

012/64 

N205 

J  12/70 

N3 

0  3/77  0 

J 12/70  OH 

0  3/77 

02 

0  6/61 

03 

OF 

“  0.0 

EFFECTIVE  FUEL 

EFFECTIVE  OXIDANT 

MIXTURE 

ENTHALPY 

HP P ( 2  ) 

IIPP{  1  ) 

HSUBO 

(KG 

-MOL  )( DEG  Kl/KG 

-0. 36660385E+02 

0.0 

-0.36660385E+02 

KG- 

ATOMS/KG 

BOP ( 1 , 2  ) 

BOP ( 1,1) 

BO (  I  ) 

C 

0 . 27577896E-02 

0.0 

0 . 27577896E 

-02 

0 

0.27577896E-02 

0.0 

0 . 27577896E 

-02 

H 

0. 27577879E-02 

0.0 

0.27577879E 

-02 

N 

0.2  757  7879F-02 

0.0 

0.27577879E 

-02 

AR 

0. 2206231 7E -01 

0.0 

0 . 220623 1 7E 

-01 

PT 

C 

0 

H  N  AR 

1 

44.084 

-115.450 

-25.642  -12.980  -18.766 

15.000 

2 

-1.293 

-31 .895 

-0.120  11.726  -15.413 

15 . 000 

3 

-2.395 

-28 . 827 

-0  043  -11.735  -15.181 

5.000 

4 

-3.2?5 

-27.055 

-8.046  -11.763  -15.068 

4  ,  000 

5 

-3.915 

-25, 7B8 

-8.060  -11.790  -14.994 

3.000 

b 

-4.463 

-24.841 

-8.0/5  -i 1.81 4  14.930 

3.000 

7 

-4.901 

-24.  105 

0.091  -11.835  14.094 

3.000 

8 

-5.260 

-23.514 

P.106  -11.054  -14.897 

3 . 000 

9 

-5.560 

-23.027 

-8.119  -11.870  14.025 

3.000 

10 

-5.815 

-22.616 

-8.132  -11.805  -14.797 

3 . 000 

1  1 

-  6.036 

-22.263 

-8.145  -11.899  14.772 

3.000 

12 

-6 . 230 

-21 . 956 

-8.15G  -11.912  14.750 

3.000 

1  3 

-6.402 

-21 . 685 

-0.168  -11.924  -14.729 

3 . 000 
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THERMODYNAMIC  EQUILIBRIUM  PROPERTIES  AT  ASSIGNED 
ENTROPY  AND  PRESSURE 


CASE 

NO. 

70 

CHEMICAL  1 

FUEL 

C 

1  .00000 

FUEL 

H 

2.00000 

FUEL 

N 

2.00000 

FUEL 

AR 

1  .00000 

MOLES  ENERGY 

CAL/MOL 
0.10000  “264  16. 162 

0.05000  0.0 

0.05000  0.0 

0.B0000  0.0 


STATE 

TEMP 

DENSITY 

DEG  K 

G/CC 

G 

298. 15 

0.0 

G 

290.15 

0.0 

G 

290.15 

0.0 

G 

298.  15 

0 . 0 

O/F-  0.0 


PERCENT  FUEL-100. 0000  EQUIVALENCE  RATIO-  2.5000 


REACTANT  DENSITY-  0.0 


THERMODYNAMIC  PROPERTIES 


P.  ATM 

T.  DEG  K 

RHO,  G/CC 

H.  CAL/G 

S.  CAL /( G ) ( K ) 

l  .0000 
300 

1 .5758-3  6, 
-119.2 
1.0280 

197.39 

138! 

.3614-2  1 
85.1 

1  . 0280 

394.77 

1732 

.0096-1  1 
143.6 

1 .0280 

592.16 

1990 

.3167-1  1 
184.9 
1.0780 

789.55 

2198 

.5891-1  1 
217.8 

1 .0280 

986.93 

2374 

.8389-1 

245.7 

1.0200 

M.  MOL  WT 
(DLV/DLP )T 

/  n  i  i/  y  m  t  <  n 

38.793 
-1 . 00094  -I 

36.513 

1  .00540  - 

36.335 

1 .00187  - 

36.304 
1.00103  - 

36.296 
1.00075  - 

36.294 
! . 00066 

CP.  CAL/1G )( K) 
GAMMA  <  S ) 

SON  VEL, M/SEC 


1.0219 

0.1660 

1.4734 

307.8 


1 . 0575 
0.  1884 
1 .4657 
678.8 


1.0151 
0.1633 
1 .5226 
776.7 


1 . 0064 
0.  1597 
1  .  5294 
834.8 


1 .0034 
0.1589 
1 .5295 
877.5 


04.32 

2528 


270.2 
1  .0280 


1381.71 

2665 


292.1 
1 .0280 


1579.10 

2790 


312.0 
1 . 02B0 


1  .0024 
0 , 1590 
1 .5274 
911.4 


1776.48 

2904 


330.4 
I  .0280 


1973.87 

3010 


347.4 
1 .0280 


36.293  36.292  36.291 

“1.00065  -1.00067  -1.00071 

1.0022  1.0024  1.0029 

0.1595  0.1602  0.1609 

1.5247  1.5217  1.5186 

939.7  963.9  985.2 


36.290  36.288 

“1.00077  -1.00084 


1  .0034 
0.1617 
1.5155 
1004.2 


1.0041 
0.  1626 
1.5125 
1021.2 


2171.26 

3108 

3.0893-1 

363.4 

1.0280 

36.286 
■  1  .  00091 
1  .0049 
0.  1634 
1  .5095 
1036.8 


MOLE  FRACTIONS 


AR 

CH20 

CH3 

CH  4 

CO 

COP 

C2H2 

C2H4 

C2HG 

H 

HCN 

HCO 

H2 

H20 

NH2 

NH3 

N2 

OH 

ADDITIONAL 


C 

C20 
H20( S) 
N20 


0.85587 

0.80556 

0. B01G3 

0.30095 

0.0 

0.00000 

0.00000 

0.00000 

0.0 

0  .  OOOOO 

0.00000 

0.00000 

0.00229 

0.00343 

0 . 00008 

0.00035 

0.03703 

0.09539 

0.09866 

0.09905 

0 . 03497 

0.00182 

0.00042 

0.00025 

0.0 

0.00000 

0.00001 

0 .00002 

0.00012 

0.00000 

0.00000 

0.00000 

0.01623 

0.00000 

0.00000 

0.00000 

0.0 

0  .  OOOOO 

0 . OOOOO 

0.00001 

0.0 

0.00004 

0.00022 

0,00043 

0.0 

0.0 

0.00000 

0.00000 

0.0 

0.04175 

0.04746 

0.04051 

0.0 

0.00166 

0.00071 

0.00056 

0.0 

0.0 

0.00000 

0  .  OOOOO 

0.0 

0.00003 

0.00003 

0.00003 

0.05349 

0.05031 

0.04997 

0.04983 

0,0 

0.0 

0 . 0 

0.00000 

PRODUCTS  WHICH  WERE  CONSIDERED  BUT  WHOSE 


CH 

CH2 

CN 

C  3 

C302 

C  4 

H20( L ) 

H202 

N 

N204 

N205 

NT 

0.80078 
0.00000 
0.00000 
0.00017 
0.09908 
0.00022 
0.00002 
0.00000 
0 . 00000 
0. 00004 
0.00057 
0.00000 
0.04876 
0.00057 
0.00000 
0.00003 
0.04975 
0.00000 


0.80073 
0.00000 
0.00000 
0.00010 
0 . 09905 
0.00021 
0.00002 
0.00000 
0.0 

0.00010 

0.00066 

0.00000 

0.04878 

0.00060 

0.00000 

0.00003 

0.04970 

0.00000 


0.80071 
0.00000 
0.00001 
0.00007 
0.09900 
0 . 00021 
0.00003 
0.00000 
0.0 

0.00018 
0.00074 
0 . 00000 
0.04B71 
0.00065 
0 . 00000 
0.00003 
0.04966 
0.00000 


0.80069 
0.00000 
0.00001 
0.00005 
0.09894 
0.00022 
0 . 00003 

o. onooo 
0.0 

0.00029 
0.00080 
0.00000 
0.04061 
0.00070 
0.00000 
0.00003 
0 . 04963 
0 . 00000 


0.80067 
0.00000 
0.00001 
0 . 00004 
0 . 09880 
0.00022 
0.00003 
0.00000 
0.0 

0;  00043 
0 . 00006 
0.00001 
0.04847 
0.00074 
0.00000 
0 . 00003 
0 .04959 
0  .  OOOOO 


0.80064 
0.00000 
0 . 00001 
0 . 00003 
0.09882 
0.00023 
0.00003 
0.00000 
0.0 

0.00060 
0.00091 
0.00001 
0.04832 
0 . 00079 
0.00000 
0.00003 
0.04957 
0.00000 


0.80060 

0.00001 

0.00001 

0.00003 

0.09875 

0.00024 

0.00003 

0.00000 

0.0 

0.00080 

0.00096 

0.00001 

0.04816 

0.00083 

0.00000 

0.00003 

0.04954 

0.00000 


0. 800S5 
0.00001 
0.00001 
0 . 00002 
0 . 09869 
0.00024 
0.00003 
0.00000 
0.0 

0.00102 
0.00100 
0.00002 
0.04798 
0.00087 
0.00001 
0 . 00003 
0.04951 
0 . 00001 


MOLE  FRACTIONS  WERE  LESS  THAN  Q.50000F-05  FOR  ALL  ASSIGNED  CONDITIONS 


CNN 

CN2 

C2 

C5 

HNCO 

HNO 

NCO 

NH 

NO 

0 

02 

03 

C2H  C2N  C2N2 

HNO?  HN03  HO 2 

N02  N03  N2II4 


PT 

C 

0 

II 

N 

AR 

1 

-6 . 556 

-21  .444 

-8.170 

-1 1.934 

-14.710 

9.000 

2 

-6.695 

-21 . 227 

-0.180 

-11.944 

-U.fi'12 

3.000 

3 

-6,821 

-21 .031 

-0.190 

-1 1.954 

-14.675 

3 . 000 

4 

-6 . 937 

-20.052 

-0.207 

-11 . 9G3 

-14,659 

3 . 000 

5 

-7.043 

-20.609 

-IJ.21G 

-11.971 

-14.645 

3.000 

6 

-7.142 

-20.530 

-0.22b 

-11.979 

- 1 4 . 6  ?l 

3.000 

7 

-7.233 

-20.399 

-8.234 

1 1 . 9CG 

-14.617 

3.000 

8 

-7.318 

-20.269 

-0.242 

-11.993 

-14.605 

3.000 

9 

-7.397 

-20  .  149 

-0.250 

-11 . 999 

14.592 

3 .000 

10 

-7.471 

-20.036 

-0.258 

-12.006 

-14.501 

3.000 

1  1 

-7.541 

-19.930 

-0.265 

12.012 

14. 570 

3.000 

12 

-  7 . GO  7 

-19.031 

8.273 

-12.017 

-14.559 

3.000 

13 

-7.669 

-19.737 

-0.200 

-  12.023 

-14.540 

.3  .  OOO 

2368.64 

3200 
.2728-1 
378.4 
l .0280 

36.283 
1  .00100 
1  .  0053 
0.1643 
1 .5066 
1051  .  J 


0 . 80049 
0 . 00001 
0 . 0000  1 
0 . 00002 
0 . 09362 
0.00025 
0 . 00004 
0 . 00009 
0 . 0 

0.0012? 
0.00105 
0 . 00002 
0.04789 
0.00091 
0. 0000 l 
0.0000  1 
0 .0  494  3 
0 , 00001 
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THERMODYNAMIC  EQUILIBRIUM  PROPERTIES  AT  ASSIGNED 


ENTROPY  AND  PRESSURE 


CASE 

NO.  70 

MOLES 

ENERGY 

STATE 

TEMP 

DENSITY 

CHEMICAL  FORMULA 

CAL/MOL 

DEG  K 

G/CC 

FUEL 

C  1.00000  0  l. 00000 

0.10000 

-264 16.152 

G 

298.15 

0.0 

FUEL 

H  Z. 00000 

0.05000 

0.0 

G 

298.15 

0.0 

FUEL 

N  2.00000 

0.05000 

0.0 

G 

298.15 

0.0 

FUEL 

AR  1.00000 

0.80000 

0.0 

G 

298.15 

0.0 

O/F5*  0.0  PERCENT  FUEL=1 00 . 0000  EQUIVALENCE  RATIO*  2.S000  REACTANT  DENSITY*  0.0 


THERMODYNAMIC  PROPERTIES 


p , 

ATM 

2566.03 

2763.42 

2960.80 

315B.  19 

3355.50 

3552.97 

3750.35 

3947.74 

4145.12 

4342.51 

4539.90 

4737.29 

4934.67 

T. 

DEG  K 

3287 

3369 

344  7 

3522 

3593 

3660 

3726 

3788 

3049 

3907 

3964 

4018 

4071 

RHO,  G/CC 

3.4515-1 

3.6261-1 

3.7960-1 

3.9640-1 

4.1281-1 

4.2892-1 

4,4477-1 

4.6036-1 

4.7573-1 

4.9087-1 

5.0581-1 

5.2056-1 

5.3513-1 

H. 

CAL/G 

392 . 7 

406.2 

419.0 

431  .4 

443.2 

454.5 

465.5 

476.0 

486.3 

496  .  1 

505 . 7 

515.1 

524  .  1 

s , 

CAL / ( G ) ( K ) 

1  .0280 

1.0280 

1 .0280 

1.0280 

1.0280 

1.0280 

1.0280 

1 .0280 

1.0280 

1 .0280 

1 .0280 

1 .0280 

1 .0280 

M, 

MOL  UT 

36.280 

36.277 

36.273 

36.269 

36.265 

36.261 

36.256 

36.251 

36.247 

36.242 

36.237 

36 . 232 

36.22/ 

( DLV/DLP  )T 

-1  .00108 

-1  .00117 

-1.00127 

-1.00137 

-1.00147 

-1 .00157 

-I .00167 

-1 .00178 

-1  .00188 

-1 .00199 

-1.00209 

-1 .00220 

-1  .00231 

( DLV/DLT  )P 

1.0067 

1 .0076 

1  .0086 

1 .0095 

1 .0105 

1.0115 

1.0125 

1 .0135 

1.0144 

1 .0154 

1.0164 

1.0173 

1.0182 

CP  , 

,  CAL/1GHK) 

0.1651 

0.1660 

0. 1669 

0.1677 

0.1605 

0. 1694 

0.1702 

0.1709 

0.1717 

0. 1724 

0.1732 

0. 1739 

0.1745 

GAMMA  IS) 

1.5038 

1.5011 

1 . 4985 

1 .4961 

1 . 4937 

1.4914 

1  .  4893 

1 .4872 

1 .4853 

1  .  4834 

1  .4816 

1  . 4799 

1 .4783 

SON  VEL. M/SEC 

1064.3 

1076.6 

1088 . 1 

1099.0 

1109.2 

1118.8 

1128.0 

1136.8 

1145,1 

1153.1 

1160.8 

1168.2 

1175.3 

MOLE  FRACTIONS 


AR 

0.90042 

0.80035 

0.80027 

0.80018 

0.80009 

0.80000 

0.79990 

0.79979 

0.79969 

0.79958 

0.79947 

0.79936 

0.79925 

CH2 

0.00000 

0 . 00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00001 

0 . 00001 

0.00001 

CH20 

0.00001 

0.00001 

0.00001 

0.00001 

0.00001 

0.00001 

0.00001 

0.00001 

0.00001 

0.00001 

0 . 0000  1 

0.00001 

0.0000  I 

CH3 

0.00001 

0.0000  1 

0.00001 

0.00002 

0.00002 

0.00002 

0.00002 

0.00002 

0.00002 

0.00002 

0.00002 

0 . 00002 

0 . 00002 

CH4 

0.00002 

0,00002 

0.00002 

0.00001 

0.00001 

0.00001 

0.00001 

0.00001 

0.00001 

0.00001 

0.00001 

0.00001 

0 . 0000  1 

CN 

0.00001 

0.00001 

0.00001 

0.00001 

0.00001 

0.00002 

0.00002 

0.00002 

0.00003 

0.00003 

0 . 00004 

0 . 00004 

0 . 00005 

CO 

0 . 09856 

0.09849 

0.09043 

0.09836 

0.09830 

0.09823 

0.09816 

0.09810 

0.09803 

0.09797 

0 . 09790 

0.09783 

0 . 09777 

C02 

0.00025 

0.00026 

0.00027 

•  0.00027 

0.00028 

0 . 00028 

0.00029 

0.00029 

0.00030 

0.00031 

0.00031 

0.00032 

0 . 00032 

C2H 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00001 

0.00001 

0.00001 

o.ooooi 

0.00001 

0.00001 

0.00001 

0 . 0000  1 

C2H2 

0.00004 

0.00004 

0.00004 

0.00004 

0.00004 

0.00004 

0.00005 

0.00005 

0.00005 

0.00005 

0.00005 

0.00005 

0.00005 

H 

0.00151 

0.00179 

0.00208 

0.00238 

0.00269 

0.00301 

0.00334 

0 . 00367 

0.00401 

0.00435 

0.00470 

0.00505 

0.00540 

HCN 

0.00109 

0.00113 

0.00117 

0.00120 

0.00124 

0.00127 

0.00130 

0.00133 

0.00136 

0.00139 

0.00142 

0.00145 

0.0014W 

HCO 

0.00002 

0.00003 

0.00003 

0.00004 

0.00004 

0.00005 

0.00005 

0.00006 

0.00006 

0.00007 

0.00008 

0.00008 

0 . 00009 

HNCO 

0.00000 

0.00000 

0.00000 

0.00000 

0.00001 

0.00001 

0.00001 

0.00001 

0.00001 

0.00001 

0.00001 

0.00001 

0.0000  1 

H2 

0.04760 

0.04740 

0.04719 

0.04697 

0.04675 

0.04653 

0.04631 

0.04608 

0.04585 

0.04562 

0.04538 

0.04515 

0.0449' 

H20 

0.00094 

0.00098 

0.00101 

0.00104 

0.00107 

0.00110 

0.00112 

0.00115 

0.00118 

0.00120 

0.00122* 

0.00124 

0.00127 

N 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0 . 00000 

0.00001 

0.00001 

0 . 0000 1 

NH 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.0000  1 

0 . 00001 

0 . 0000  1 

NH2 

0.00001 

0 . 00001 

0.00001 

0.00001 

0.00002 

0.00002 

0.00002 

0.00002 

0 . 00003 

0.00003 

0 . 00003 

0.00003 

0 . 00004 

NH3 

0.00003 

0.00004 

0.00004 

0.00004 

0.00004 

0.00004 

0.00004 

0.00004 

0.00004 

0.00004 

0.00004 

0.00004 

0.0000  1 

NO 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00001 

0.00001 

0.00001 

0.00001 

0.00001 

0.00001 

0.00002 

0.00002 

N2 

0.04945 

0.04943 

0 . 04940 

0.04937 

0.04935 

0.04932 

0.04929 

0.04926 

0.04924 

0.04921 

0.04918 

0.04915 

0 .04913 

0 

0.00000 

0,00000 

0  .00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00001 

0.0000 1 

OH 

0.00001 

0.00002 

0.00002 

0.00002 

0.00003 

0.00004 

0.00004 

0.00005 

0.00006 

0.00007 

0.00007 

0 . 00008 

0.0000' 

ADDITIONAL 

PRODUCTS  WHICH  WERE  CONSIDERED 

BUT  WHOSE 

MOLE  FRACTIONS  WERE  LESS  THAN  0.50000E-05  FOR  ALL  ASSIGNED  CONDITIONS 

C 

CH 

CNN 

CN2 

C2 

C2H4 

C2H6 

C2N 

C2N2 

C20 

C  3 

C302 

C  4 

C5 

HMO 

HNOZ 

HN03 

H02 

H20< S  > 

H20( L  ) 

H202 

NCO 

NO  2 

N03 

N2IM 

N20 

N204 

N205 

N3 

02 

03 
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